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A new critical survey is presented of all half-life, decay-energy and branching-ratio measurements 
related to 20 superallowed + — > + (3 decays. Compared with our last review, there are numerous 
improvements: First, we have added 27 recently published measurements and eliminated 9 refer- 
ences, either because they have been superceded by much more precise modern results or because 
there are now reasons to consider them fatally flawed; of particular importance, the new data in- 
clude a number of high-precision Penning-trap measurements of decay energies. Second, we have 
used the recently improved isospin symmetry-breaking corrections, which were motivated by these 
new Penning-trap results. Third, our calculation of the statistical rate function / now accounts 
for possible excitation in the daughter atom, a small effect but one which merits inclusion at the 
present level of experimental precision. Finally, we have re-examined the systematic uncertainty as- 
sociated with the isospin symmetry-breaking corrections by evaluating the radial-overlap correction 
using Hartree-Fock radial wave functions and comparing the results with our earlier calculations, 
which used Saxon- Woods wave functions; the provision for systematic uncertainty has been changed 
as a consequence. The new "corrected" Tt values are impressively constant and their average, 
when combined with the muon liftime, yields the up-down quark-mixing element of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, V ud = 0.97425 ± 0.00022. The unitarity test on the top row of 
the matrix becomes |Vud| 2 + Vus] 2 + Vuf,| 2 = 0.99995 ± 0.00061. Both V u d and the unitarity sum 
^ | have significantly reduced uncertainties compared with our previous survey, although the new value 

D , of V u d is statistically consistent with the old one. From these data we also set limits on the possible 

existence of scalar interactions, right-hand currents and extra Z bosons. Finally, we discuss the 
priorities for future theoretical and experimental work with the goal of making the CKM unitarity 
^ , test even more definitive. 
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PACS numbers: 23.40.Bw, 12.15.Hh, 12.60.-i 



INTRODUCTION 



\ Precise measurements of the beta decay between nuclear analog states of spin, J~* = + , and isospin, T = 1, provide 
demanding and fundamental tests of the properties of the electroweak interaction. Collectively, these transitions can 
sensitively probe the conservation of the vector weak current, set tight limits on the presence of scalar or right- 
hand currents and, by providing the most precise value for V uc i, the up-down quark- mixing element of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, they can contribute to the most demanding available test of the unitarity of that 
• i-H , matrix, a property which is fundamental to the electroweak standard model. 

We have published five previous surveys of + superallowed transitions [l|, 0, [H, EL IE] , the first having appeared 
35 years ago and the most recent, four years ago. In each, we published a complete survey of all relevant nuclear data 
that pertained to these superallowed transitions and used the results to set limits on the weak-interaction parameters 
that were important at the time. A particularly noteable outcome of our analysis four years ago Q was that the sum 
of squares of the top-row elements of the CKM matrix - the test of CKM unitarity - remained ambiguous, with the 
possibility of a significant shortfall in the unitarity sum. 

Since our last survey closed in November 2004, there has been a great deal of activity in this field prompted at 
least in part by the tantalizing possibility that new physics could be revealed by a failure in CKM unitarity. New 
measurements relating to + — > + superallowed transitions have appeared in 27 publications, an addition of 20% to 
the papers accumulated up to 2004. Many of these measurements were of unprecedented precision so they did not 
merely add more of the same: they palpably improved the results, in some cases by tightening their error bars and, in 
others, by changing their central values. Penning-trap measurements of decay energies, which only became possible 
after 2004, have been especially effective in this regard. 

In addition to new measurements, there have also been important improvements to the small theoretical corrections 
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that must be applied to the data in order to extract V u d and the values of other weak-interaction parameters. In the 
past four years, the radiative Q and isospin symmetry-breaking corrections 0] have both been subjected to major 
re-evaluations, which have undoubtedly improved their values and, in the former case, has reduced the uncertainty 
by a factor of two. 

In parallel with these developments, there has also been considerable activity in the determination of V US} the other 
matrix element that plays a role in the top-row unitarity test of the CKM matrix. (The third element in the top row, 
V u b, is very small and contributes a negligible 0.001% to the unitarity sum.) As with the work related to V u d, this 
activity has encompassed new experiments - precise measurements of kaon branching ratios - as well as improved 
theoretical corrections. However, in contrast with VyA, not only the uncertainty of V us but also its central value have 
been considerably changed by this recent work (see Q for an up-to-date overview of V us ). 

Overall, the recent improvements have been numerous enough and their impact on the unitarity test significant 
enough that this is an opportune time to produce a new and updated survey of the nuclear data used to establish V uc i- 
We incorporate data on a total of 20 superallowed transitions and have continued the practice we began in 1984 Q 
of updating all original data to take account of the most modern calibration standards. In addition to including the 
improved correction terms already referred to, we have also upgraded our calculation of the statistical rate function / 
to include provision for excitation of the daughter atom, and we have included a more extensive treatment of possible 
systematic uncertainties associated with the isospin symmetry-breaking corrections. 

Superallowed + ^ + (3 decay between T — 1 analog states depends uniquely on the vector part of the weak 
interaction and, according to the conserved vector current (CVC) hypothesis, its experimental ft value should be 
related to the vector coupling constant, a fundamental constant which is the same for all such transitions. In practice, 
the expression for ft includes several small (~1%) correction terms. It is convenient to combine some of these terms 
with the ft value and define a "corrected" Tt value. Thus, we write @ 

Tt = ft(l + ^)(1 + S N8 - 5c) = 2G z* + , (1) 

where K/(hc) 6 = 27r 3 ^ln2/(m e c 2 ) 5 = 8120.2787(11) x 10~ 10 GeV~ 4 s, G v is the vector coupling constant for semi- 
lcptonic weak interactions, 8c is the isospin-symmetry-breaking correction and is the transition-independent part 
of the radiative correction. The terms 8' R and 8ns comprise the transition-dependent part of the radiative correction, 
the former being a function only of the electron's energy and the Z of the daughter nucleus, while the latter, like 8c, 
depends in its evaluation on the details of nuclear structure. From this equation, it can be seen that each measured 
transition establishes an individual value for G v and, if the CVC assertion is correct that G v is not renormalized in 
the nuclear medium, all such values - and all the Tt values themselves - should be identical within uncertainties, 
regardless of the specific nuclei involved. 

Our procedure in this paper is to examine all experimental data related to 20 superallowed transitions, comprising all 
those that have been well studied, together with others that are now coming under scrutiny after becoming accessible 
to precision measurement. The methods used in data evaluation are presented in Sec. [TTJ The calculations and 
corrections required to extract Tt values from these data are described and applied in Sec. IIII1 in the same section, 
we use the resulting Tt values to test CVC. Finally, in Sec. IIVI we explore the impact of these results on a number 
of weak-interaction issues: CKM unitarity as well as the possible existence of scalar interactions, right-hand currents 
and extra Z bosons. This is much the same pattern as we followed in our last review [H so we will not describe the 
formalism again in detail, referring the reader instead to that earlier work. 

II. EXPERIMENTAL DATA 

The /i-value that characterizes any /3 transition depends on three measured quantities: the total transition energy, 
QeCi the half-life, ti/2, of the parent state and the branching ratio, R, for the particular transition of interest. The 
Q.EC- v alue is required to determine the statistical rate function, /, while the half-life and branching ratio combine to 
yield the partial half-life, t. In Tables Hl lVIII we present the measured values of these three quantities and supporting 
information for a total of twenty superallowed transitions. 

A. Evaluation principles 

In our treatment of the data, we considered all measurements formally published before September 2008 and those 
we knew to be in an advanced state of preparation for publication by that date. We scrutinized all the original 
experimental reports in detail. Where necessary and possible, we used the information provided there to correct 
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FIG. 1: Differences between individual measurements and the averages of all measurements for the seven parent nuclei studied 
by Vonach et al. [Vo77]. The filled squares are the results of the ( 3 He, t) measurements of Vonach et al.; the open squares are 
recent Penning-trap results [Sa05, Er06b, Er08, Ge08]. For each parent nucleus, the grey band about the zero line represents 
the uncertainty of the average for that case. Note that all the averages include the results of Vonach et al., the Penning-trap 
results and any other relevant measurements appearing in Table U 



the results for calibration data that have improved since the measurement was made. If corrections were evidently 
required but insufficient information was provided to make them, the results were rejected. Of the surviving results, 
only those with (updated) uncertainties that are within a factor of ten of the most precise measurement for each 
quantity were retained for averaging in the tables. Each datum appearing in the tables is attributed to its original 
journal reference via an alphanumeric code comprising the initial two letters of the first author's n ame and the two last 
digits of the publication date. These codes are correlated with the actual reference numbers, (9j- [l73j |. in Table IVTTT1 
The statistical procedures we have followed in analyzing the tabulated data are based on those used by the Particle 



Data Group in their periodic reviews of particle properties (e.g. Ref. 174 1) and adopted by us in our previous surveys. 
We gave a detailed description of those procedures in our 2004 survey [5] so will not repeat it here. 

Our evaluation principles and associated statistical procedures constitute a very conservative approach to the data. 
Unless there is a clearly identifiable reason to reject a result, we include it in our data base even if it deviates 
significantly from other measurements of the same quantity, the consequent non-statistical spread in results being 
reflected in an increased uncertainty assigned to the average. Wherever this occurs, the factor by which the uncertainty 
has been increased is listed in the "scale" column of a table. There are a few examples, though, of a single publication 
that includes a number of measurements - a set of half- lives or Qec values for example - most or all of which 
deviate substantially from other accepted measurements of the same quantities. In such cases, we consider that some 
systematic problem has been revealed, and exclude all the results from that publication. If any measurement with an 
acceptable uncertainty is nevertheless excluded from our data base, the reason for its exclusion is listed in Table IVIII 

One particlarly significant, longstanding reference had to be excluded for the first time from this sur vey. Our 
decision to do so deserves a more detailed explanation. In 1977, Vonach et al. published in a single paper [162j | the 
Qec values for seven superallowed emitters ( 14 0, 26 Al m , 34 C1, 42 Sc, 46 V, 50 Mn and 54 Co), which they had determined 
from the Q values for ( 3 He,i) reactions on their stable daughters. They had used a "precision time-of- flight measuring 
system" with the Q3D spectrograph of the Munich Tandem Laboratory to produce uncertainties of 0.4-0.6 keV. For 
the time, these were very precise results and consequently they had a major impact on the superallowed data base 
for the following three decades. 

The first indication that the Vonach et a l. results might have a problem came with the first Penning-trap mea- 
surement of a superallowed Qec value [l49j |. The new measurement for 46 V quoted 0.4-keV uncertainty and differed 
from the old result by 2.4 keV, four of Vonach's claimed standard deviations. Within a year, a second Penning-trap 
measurement [6l| had confirmed the new 46 V result and had also found that the 42 Sc Qec value differed from the 
Vonach result by six times the latter's quoted uncertainty. Two years later, another Penning-trap measurement [r32| 
found the 50 Mn and 54 Co Qec values also differed from Vonach's results by a similar amount. These most-recent 
Penning-trap results quoted 0.1-keV uncertainties. A current overview of the situation for all seven of the superallowed 
transitions measured by Vonach appears in Fig. [TJ where each Vonach result is compared with the equivalent result (s) 
from a Penning trap, and both are compared to the average of all results for the same transition. With only two 
transitions, those from 26 Al m and 34 C1, showing agreement and the four cases already mentioned displaying serious 
disagreement, w e be lieved that the best approach was to eliminate all the results published in the original Vonach 
et al. reference jl62l |. This conclusion has been further supported by a very recent ( 3 He,i) measurement of the 46 V 
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FIG. 2: For the "traditional nine" transitions, we plot the differences between individual measurements and average values; 
the results of Vonach et al. [Vo77] have been removed from the averages. The open squares are the results of Penning-trap 
measurements [Sa05, Er06b, Er08, Ge08]; the filled squares are from combined (p, 7) and (n, 7) measurements [De69 and the 
references listed in footnotes 6, 7 and 8 of Table [J; the filled triangles are from (p,n) threshold measurements [Ba84, Ba98, 
Wh77, To03, Br94, Ba77c, Ja78, Ha98]; and the filled diamond is the new ( 3 He,i) measurement for 46 V [Fa09]. For each parent 
nucleus, the grey band about the zero line represents the uncertainty of the average for that case. 



Qec value [631 ] made with much of the same experimental equipment originally used by Vonach et al. thirty years 
ago. The new result disagrees with the old measurement and confirms the new Penning-trap values. 

Before Penning traps could be applied to these measurements, all superallowed Qec values were determined via 
nuclear reactions. In addition to those employing ( 3 He,i) reactions, two other types of experiment led to rather 
precise results: the measurement of (p, n) thresholds and the combined measurements of (p, 7) and (n, 7) Q~values 
on the same target, one reaction leading to the superallowed parent and the other to the daughter. We are now in 
a position to compare the different types of m easurement to examine whether there are any systematic differences 
among them. A careful study of this issue [l75j . restricted to the region around A=26, was undertaken several years 
ago and found no evidence of any systematic differences. We can now confirm this conclusion over a wider mass 
range with the help of Fig. [5] In that figure we consider nine superallowed transitions, which we will refer to as the 
"traditional nine" cases. They are the only superallowed transitions that populate a stable daughter nucleus and, 
for obvious reasons, were the only ones whose Q values could be measured to high precision in the pre-trap era. 
There are no systematic deviations app arent in the figure, leading us to conclude that, whatever problems plagued 
the measurements of Vonach et al. (1621 ] , they were associated with that particular experiment and were not endemic 
to a whole class of experiments. Of course, this conclusion could be strengthened by new Penning-trap data for 10 C, 
14 0, 34 C1 and 38 K m 



B. Data Tables 



TABLE I: Measured results from which the decay transition energies, Qec-, have been derived for superallowed /3-dccays. The lines 
giving the average superallowed Qec values themselves are in bold print. (Sec Tabic |VI III for the correlation between the alphanumeric 
reference code used in this table and the actual reference numbers.) 



Parent/Daughter Property 1 Measured Energies used to determine Qec (kcV) Average value 

nuclei 1 2 3 Energy (kcV) scale 



Qec(ss) 


3647.84 ± 0.34 


[Ba84] 


E m (d0+) 


1740.15 ± 0.17 


[AJ88] 


Qec(so) 






QEc(gs) 


5143.35 ± 0.60 


[Bu61] 




5143.43 ± 0.37 [Wh77] 


2312.798 ± 0.011 [Aj91] 


Qec(scl) 






ME(p) 


5316.8 ± 1.5 


[Ma94] 


ME(d) 


873.31 ± 0.94 


[Bo64] 




877.2 ± 3.0 


[Sc73] 


QEc(gs) 


4438 ± 9 


[Fr63] 


E a; (d0+) 


1041.55 ± 0.08 


[Ti95] 









3647.95 ± 0.12 
1740.07 ± 0.02 



5145.09 ± 0.46 [Ba62] 
5144.34 ± 0.17 [To03] 



5317.63 ± 0.36 [B104b] 
875.5 ± 2.2 [Ho64] 
873.96 ± 0.61 [Ro75] 



5145.57 ± 0.48 [Ro70] 



876.5 ± 2.8 [Pr67] 



3647.94 ± 0.11 
1740.07 ± 0.02 
1907.87 ± 0.11 



5144.33 ± 0.29 
2312.798 ± 0.011 
2831.24 ± 0.23 ! 

5317.58 ± 0.35 

874.02 ± 0.48 

4443.54 ± 0.60 

1041.55 ± 0.08 
3401.99 ± 0.60 



1.0 
1.0 



1.0 



1.0 
1.0 
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TABLE I (continued) 



2 Mg 



"Si 



3 Ca 



T z = 0: 



Parent/Daughter Property 1 


Measured Energy (keV) 




Average value 


nuclei 


1 2 


3 


Energy (keV) scale 



"Al 



ME{p) 
ME(d) 

QEc(gs) 

E x (d0+) 
Qec(scl) 

ME(p) -7145.4 ± 3.0 [Ha74c 

ME(dO + ) -11981.99 ± 0.26 5 
Qsc(sa) 4850 ± 13 [Fr63] 



-401.2 ± 3.0 
-5184.3 ± 1.5 
-5183.2 ± 1.0 
4781.64 ± 0.28 

657.00 ± 0.14 



[Ha74c] 

[Wo68] 

[Gi72] 

[Mu04] 

[En98] 



ME{p) 

ME(d) 
QEc(gs) 
E a: (d0 + ) 
Qec(scl) 



-14060 ± 15 
-14063.4 ± 3.0 
-20203 ± 3 

677.29 ± 0.07 
5437 ± 17 



[Mi67] 

[Ha74c] 

[Ha67] 

En98] 
Rr63] 



CI 



ME{p) -18380.2 ± 3.0 [Ha74c] 
ME(d) -24440.15 ± 0.26 5 
Qsc(sa) 



ME(p) 
ME(dO + ) 

ME(p) 
ME(d) 
Qec(so) 



-22058.53 ± 0.28 [Ri07] 
-28670.20 ± 0.32 5 



-25121 ± 6 [Mi67] 
-32121.12 ± 0.29 5 



-400.4 ± 1.3 4 -400.5 ± 1.0 [Pa05] 

-5182.5 ± 0.5 [Bc68] -5181.3 ± 1.7 [An70] 

-5181.56 ± 0.16 [Mu04] -5181.08 ± 0.30 [Sa04] 
4781.40 ± 0.67 [Sa04] 



-7139.5 ± 1.0 [Pa05] 

-14054 ± 25 [Mc67] -14068 ± 30 [Ha68] 
-20200.58 ± 0.40 [Re85] 

-18378.4 ± 3.5 [HoOl] -18377.10 ± 0.41 [He02] 
-22058.01 ± 0.65 [Ge07] 

-25086 ± 30 [Ha68] -25124 ± 13 [Zi72] 



-400.5 ± 0.8 

-5181.58 ± 0.19 
4781.55 ± 0.25 
657.00 ± 0.14 
4124.55 ± 0.28 

-7140.1 ± 1.8 
-11981.99 ± 0.26 
4842.0 ± 1.8 



-14063.1 ± 2.9 
-20200.62 ± 0.40 
6137.5 ± 2.9 
677.29 ± 0.07 
5459.5 ± 3.9 

-18377.17 ± 0.40 
-24440.15 ± 0.26 
6062.98 ± 0.48 

-22058.45 ± 0.26 
-28670.20 ± 0.32 
6611.75 ± 0.41 

-25120.7 ± 5.3 
-32121.12 ± 0.29 
7000.5 ± 5.4 



1.0 



1.7 
1.0 



1.9 
1.0 



1.0 
1.0 



1.3 

1.0 



1.0 



1.0 



26 M m 


26 Mg 


QEc(gs) 

E*(p0+) 
Qsc(sa) 


4004.79 ± 0.55 
228.305 ± 0.013 
4232.19 ± 0.12 


[Dc69] 
: [En98] 
[Br94] 


4004.41 ± 0.10 6 
4232.83 ± 0.13 [Er06b] 


4004.40 ± 0.22 [Ge08] 


4004.42 ± 0.09 
228.305 ± 0.013 
4232.66 ± 0.06 3 


1.0 
1.0 


34 CI 


34g 


Qsc(sa) 


5490.3 ± 1.9 
5491.65 ± 0.26 


[Ry73a] 

7 


5491.6 ± 2.3 [Ha74d] 


5491.71 ± 0.54 [Ba77c] 


5491.64 ± 0.23 


1.0 


33 K m 


38 Ar 


G>Ec(gs) 
Qsc(sa) 


5914.82 ± 0.61 
130.4 ± 0.3 
6044.6 ± 1.5 


[Ja78] 
[En98] 
[Bu79] 


6044.38 ± 0.12 [Ha98] 




5914.82 ± 0.61 
130.4 ± 0.3 
6044.40 ± 0.11 


1.0 


42 Sc 


42 Ca 


Qsc(sa) 


6425.84 ± 0.17 


8 


6426.13 ± 0.21 [Er06b] 




6426.28 ± 0.30 3 


3.0 


46y 


46 Ti 


Qsc(sa) 


7053.3 ± 1.8 
7052.11 ± 0.27 


[Sq76] 
[Fa09] 


7052.90 ± 0.40 [Sa05] 


7052.72 ± 0.31 [Er06b] 


7052.49 ± 0.16 


1.3 


50 Mn 


50 Cr 


Qec{so) 


7634.48 ± 0.07 


[Er08] 






7634.45 ± 0.07 3 


1.0 


54 Co 


54 Fo 




8244.54 ± 0.10 


[Er08] 






8244.37 ± 0.28 3 


3.4 


62 Ga 


62 Zn 


Qec(so) 


9181.07 ± 0.54 


[Er06a] 






9181.07 ± 0.54 




66 As 


66 Gc 


ME(p) 
ME(d) 
Qec(so) 


-52018 ± 30 
-61607.0 ± 2.4 
9550 ± 50 


[Sc07] 
[Sc07] 
[Da80] 






52018 ± 30 
-61607.0 ± 2.4 
9579 ± 26 


1.0 


70 Br 


70 So 


Qec(sci) 


9970 ± 170 


[Da80] 






9970 ± 170 




74 Rb 


74 Kr 


ME{p) 
ME(d) 
Qec(scl) 


-51905 ± 18 
-62332.0 ± 2.1 


[He02] 
[Ro06] 


-51915.2 ± 4.0 [Ke07] 




-51914.7 ± 3.9 
-62332.0 ± 2.1 
10417.3 ± 4.4 


1.0 



1 Abbreviations used in this column are as follows: "gs", transition between ground states; "sa", supcrallowcd transition; "p", parent; 
"<F , daughter; "ME", mass excess; U E X (0 + )" , excitation energy of the + (analog) state. Thus, for example, "Q Bc( sa )" signifies the 
Qs^-value for the superallowed transition, "M_E(d)", the mass excess of the daughter nucleus; and "ME(d0 + ), the mass excess of the 
daughter's + state. 

2 Result based on references [Ba88] and [Ba89]. 

3 Average result includes the results of Qec pairs; see Table HT1 

4 Result based on references [Bi03], [Se05] and [Je07]. 

5 Result obtained from data elsewhere in this table. 

6 Result based on references [Is80], [A182], [Hu82], [Be85], [Pr90], [Ki91] and [Wa92]. 

7 Result based on references [Wa83], [Ra83] and [Li94]. 

8 Result based on references [Zi87] and [Ki89]. 
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TABLE II: Qfic-value differences for superallowed /3-decay branches. These data are also used as input to determine some of 
the average Q^c-values listed in Table [T] (See Table [VTTT1 for the correlation between the alphabetical reference code used in 
this table and the actual reference numbers.) 



Parent 


Parent 


QeC2 — 


Qeci (keV) 




nucleus 1 


nucleus 2 


measurement 




average" 


14q 
26 A 

42 Sc 
42 Sc 
50 Mn 


26 js^y 71 
42 Sc 
50 Mn 
54 Co 
54 Co 


1401.68 ± 0.13 [Ko87] 

2193.5 ± 0.2 [Ko87] 

1207.6 ± 2.3 [Ha74d] 
1817.2 ± 0.2 [Ko87] 
610.09 ± 0.17 




1401.43 ± 0.24 
2193.62 ± 0.30 
1208.17 ± 0.30 
1818.09 ± 0.41 
609.92 ± 0.29 


"Average 


values include the results of direct Qec 


-value measurements: see Table [T] 







The Qec-value data appear in Tables HI and HT1 For the "traditional nine" superallowed decays - those of 10 C, 14 0, 
26 Al m , 34 C1, 38 K m , 42 Sc, 46 V, 50 Mn and 54 Co - with stable daughter nuclei, their Q EC values were all determined 
in the past by direct reaction measurements of that property. More recently, a growing number of Penning-trap 
measurements, also extending to nuclei outside of the traditional nine, determine the parent and daughter masses in 
a single experiment, thus effectively measuring the Qec value directly. Measurements of both types are identified 
in column 3 of Table [J by "Qbc^a)" and each individual result is itemized with its appropriate reference in the 
next three columns. The weighted average of all measurements for a particular decay appears in column 7, with the 
corresponding scale factor (see Sec. Ill A[) in column 8. A few of these cases, like 34 C1 and all the cases from 42 Sc 
to 62 Ga, have no further complications. There are other cases, however, in which Q^c-value differences have been 



TABLE III: Half-lives, ti/2, of superallowed /3-emitters. (See Table IVTnl for the correlation between the alphabetical reference 
code used in this table and the actual reference numbers.) 



Parent 



Measured half-lives, ii/2 ( ms ) 



Average value 



nucleus 






1 






2 




3 








4 


*l/2 


(ms) 


scale 


= -1: 


































io c 


19280 


± 


20 


[Az74] 


19295 ± 


15 


[Ba90] 


19310 ± 4 


[Ia08] 








19308.0 


± 


3.8 


1.0 


14 Q 


70480 


± 


150 


[A172] 


70588 ± 


28 


[C173] 


70430 ± 180 [Az74] 


70684 


± 


77 [Be78] 












70613 


± 


25 


[Wi78] 


70560 ± 


49 


[GaOl] 


70641 ± 20 


[Ba04] 


70696 


± 


52 [Bu06] 


70620 


± 


15 


1.2 


18 Ne 


1669 


± 


4 


[A175] 


1687 ± 


9 


[Ha75] 


1665.6 ± 1.9 


[Gr07] 








1667.0 


± 


1.7 


1.0 


22 Mg 


3857 


± 


9 


[Ha75] 


3875.5 ± 


1.2 


[Ha03] 












3875.2 


± 


2.4 


2.0 


26 Si 


2210 


± 


21 


[Ha75] 


2240 ± 


10 


[Wi80] 


2228.3 ± 2.7 


[Ma08] 








2228.8 


± 


2.0 


1.1 


30g 


1180 


± 


40 


[Ba67] 


1220 ± 


30 


[Mo71] 


1178.3 ± 4.8 


[Wi80] 








1179.4 


± 


4.7 


1.0 


34 Ar 


844.5 


± 


3.4 


[Ha74a] 


843.8 ± 


0.4 


[Ia06] 












843.8 


± 


0.4 


1.0 


3s Ca 


470 


± 


20 


[Ka68] 


439 ± 


12 


[Ga69] 


450 ± 70 


[Zi72] 


430 


± 


12 [W180] 


440.0 


± 


7.8 


1.2 


42 T; 


200 


± 


20 


[Ni69] 


202 ± 


5 


[Ga69] 


173 ± 14 


[A169] 








198.8 


± 


6.3 


1.4 


T z = 0: 
26 A1 ™ 


6346 


± 


5 


[Fr69a] 


6346 ± 


5 


[Az75] 


6339.5 ± 4.5 


[A177] 


6346.2 


± 


2.6 [Ko83] 


6345.0 


± 


1.9 


1.0 




6345 


± 


14 


[Sc05] 


























34 C1 


1526 


± 


2 


[Ry73a] 


1525.2 ± 


1.1 


[Wi76] 


1527.7 ± 2.2 


[Ko83] 


1526.8 


± 


0.5 [Ia06] 


1526.55 


± 


0.44 


1.0 


3s K m 


925.6 


± 


0.7 


[Sq75] 


922.3 ± 


1.1 


[Wi76] 


921.71 ± 0.65 


[Wi78] 


924.15 


± 


0.31 [Ko83] 












924.4 


± 


0.6 


[BaOO] 


924.46 ± 


0.14 [Ba08] 












924.33 


± 


0.27 


2.3 


42 Sc 


680.98 


± 


0.62 


[Wi76] 


680.67 ± 


0.28 [Ko97a] 












680.72 


± 


0.26 


1.0 


46y 


422.47 


± 


0.39 


[A177] 


422.28 ± 


0.23 [Ba77a] 


422.57 ± 0.13 


[Ko97a] 








422.50 


± 


0.11 


1.0 


50 Mn 


284.0 


± 


0.4 


[Ha74b] 


282.8 ± 


0.3 


[Fr75] 


282.72 ± 0.26 


[Wi76] 


283.29 


± 


0.08 [Ko97a] 












283.10 


± 


0.14 


[Ba06] 


















283.21 


± 


0.11 


1.7 


54 Co 


193.4 


± 


0.4 


[Ha74b] 


193.0 ± 


0.3 


[Ho74] 


193.28 ± 0.18 


[A177] 


193.28 


± 


0.07 [Ko97a] 


193.271 


± 


0.063 


1.0 


62 Ga 


115.84 


± 


0.25 


[Hy03] 


116.19 ± 


0.04 [B104a] 


116.09 ± 0.17 


[Ca05] 


116.01 


± 


0.19 [Hy05] 












116.100 


± 


0.025 [Gr08] 


















116.121 


± 


0.040 


1.9 


66 As 


95.78 


± 


0.39 


[A178] 


95.77 ± 


0.28 [Bu88] 


97 ± 2 


[Ji02] 








95.79 


± 


0.23 


1.0 


™Br 


80.2 


± 


0.8 


[A178] 


78.54 ± 


0.59 [Bu88] 












79.12 


± 


0.79 


1.7 


74 Rb 


64.90 


± 


0.09 


[OiOl] 


64.761 ± 


0.031 [BaOl] 












64.776 


± 


0.043 


1.5 
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TABLE IV: Measured results from which the branching ratios, R, have been derived for superallowed /3-transitions. The lines 
giving the average superallowed branching ratios themselves are in bold print. ( See Table ["VTTT1 for the correlation between the 
alphabetical reference code used in this table and the actual reference numbers.) 



Parent/Dau 


ghtcr 


Daughter state 


Measured Branchin 


g Ratio, R (%) 


Average value 




nuclei 




E x (MeV) 


1 


2 


R (%) 


scale 


T z = -I: 














io c 




2.16 
1.74 


+o.ooo8 [Go72] 

1.468 ± 0.014 [Ro72] 
1.465 ± 0.009 [Kr91] 
1.4665 ± 0.0038 [Fu99] 


1.473 ± 0.007 [Na91] 
1.4625 ± 0.0025 [Sa95] 


n +0.0008 
U -0 

1.4646 ± 0.0019 


1.0 


14 Q 


14 N 


gs 

3.95 
2.31 


0.68 ± 0.10 [Sh55,To05] 
0.54 ± 0.02 [Si66,T°05] 
0.062 ± 0.007 [Ka69] 
0.053 ± 0.002 [He81] 


0.74 ± 0.05 [Fr63,To05]] 
0.058 ± 0.004 [Wi80] 


0.571 ± 0.068 

0.0545 ± 0.0019 
99.374 ± 0.068 


3.7 
1.1 


ls Nc 


18 p 


1.04 


9 ± 3 [Fr63] 


7.70 ± 0.21° [Ha75] 


7.70 ± 0.21 


1.0 


22 Mg 


22 Na 


0.66 


54.0 ± 1.1 [Ha75] 


53.15 ± 0.12 [Ha03] 


53.16 ± 0.12 


1.0 


Si 


26 a i 
Al 


1.06 
0.23 


21.8 ± 0.8 [Ha75] 


21.21 ± 0.64 [Ma08] 


21.4 ± 0.5 

i D.4U ± U.o ( 


1.0 


30 g 


30p 


0.68 


20 ± 1 [Fr63] 




20 ± 1 
77.4 ± 1.0" 




34 Ar 


34 C1 


0.67 
gs 


2.49 ± 0.10 [Ha74a] 




2.49 ± 0.10 
94.45 ± 0.25° 




42 T; 


42 Sc 


0.61 
gs 


56 ± 14 [A169] 




56 ± 14 
43 ± 14° 




T z = 0: 














26 A1 ™ 


26 Mg 


gs 


>99.997 [Ki91] 




lOO.OOO+o.^ 




34 C1 


34g 


gs 


>99.988 [Dr75] 




ioo.ooot° 012 




38 K m 


38 Ar 


3.38 
gs 


<0.0019 [Ha94] 
u.uoou nr. u. uu^to [Ijcuoj 


<0.0008 [Le08] 


n +o.ooos 

u.uoou nr. u.uu^to 
99 9670 + 0043 

-0.0044 




42 Sc 


42 Ca 


1.84 


0.0063 ± 0.0026 [In77] 


0.0022 ± 0.0017 [De78] 












0.0103 ± 0.0031 [Sa80] 


0.0070 ± 0.0012 [Da85] 


0.0059 ± 0.0014 


1.6 






gs 






99.9941 ± 0.0014 




46y 


46 Ti 


2.61 
4.32 
EGT 4 
gs 


0.0039 ± 0.0004 [Ha94] 
0.0113 ± 0.0012 [Ha94] 
<0.004 




0.0039 ± 0.0004 
0.0113 ± 0.0012 

+ 0.004 

99.9848±° ;° °f 2 




Mn 


Cr 


3.63 
3.85 
5.00 
gs 


0.057 ± 0.003 [Ha94] 
<0.0003 [Ha94] 
0.0007 ± 0.0001 [Ha94] 




0.057 ± 0.003 

n +0.0003 

u -o 

0.0007 ± 0.0001 
99.9423 ± 0.0030 




54 Co 


54 Fc 


2.56 

EGT 4 

gs 


0.0045 ± 0.0006 [Ha94] 
<0.03 




0.0045 ± 0.0006 
n +0.03 

u -0 

ya - y9£>& -0.0300 




62 Ga 


62 Zn 


EGT 4 
gs 


0.142 ± 0.008 [Fi08] 


0.107 ± 0.024 [Bc08] 


0.139 ± 0.011 
99.862 ± 0.011 


1.4 


74 Rb 


74 Kr 


EGT 4 
gs 


0.50 ± 0.10 [Pi03] 




0.50 ± 0.10 
99.50 ± 0.10 





"Result also incorporates data from Table IVl 

"designates total Gamow- Teller transitions to levels not explicitly listed; values were derived with the help of calculations in [Ha02]. 



measured in addition to the individual Qec~ values. These measurements are presented in Table [TTJ They have been 
dealt with in combination with the direct Q^c-value measurements, as described in Ref. with the final average 
Qec value appearing in column 7 of Table U and the average difference in column 4 of Table HT1 Both are flagged with 
footnotes to indicate the interconnection. 

There are two cases, 26 Al m and 38 K m , in which the superallowed decay originates from an isomeric state. For both, 
there are Qsc-value measurements that correspond to the ground state as well as to the isomer. Obviously, these two 
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TABLE V: Relative intensities of /3-delayed 7-rays in the superallowed /3-decay daughters. These data are used to determine 
some of the branching ratios presented in Table HVl (See Table IVTnl for the correlation between the alphabetical reference code 
used in this table and the actual reference numbers.) 



Parent/Daughter 
nuclei 


daughter 
ratios" 






Measured 7-ray 


Ratio 








Average value 








1 






2 




Ratio 


scale 


18 Nc 18 F 


766o/7l042 


0.021 


± 


0.003 [Ha75] 


0.0169 


± 


0.0004 


[He82] 










0.0172 


± 


0.0005 [Ad83] 










0.0171 ± 0.0003 


1.0 


26 Si 26 A1 


71622 /7829 


0.149 


± 


0.016 [Mo71] 


0.134 


± 


0.005 


[Ha75] 










0.1245 


± 


0.0023 [Wi80] 


0.1301 


± 


0.0062 


[Ma08] 


0.1269 ± 0.0026 


1.3 




71655/7829 


0.00145 


± 


0.00032 [Wi80] 










0.0015 ± 0.0003 






71843/7829 


0.013 


± 


0.003 [Mo71] 


0.016 


± 


0.003 


[Ha75] 










0.01179 


± 


0.00027 [Wi80] 










0.0118 ± 0.0003 


1.0 




72512/7829 


0.00282 


± 


0.00010 [Wi80] 










0.0028 ± 0.0001 






7total/7829 
















0.1430 ± 0.0026 




30g 30 p 


7709/7677 


0.006 


± 


0.003 [Mo71] 


0.0037 


± 


0.0009 


[Wi80] 


0.0039 ± 0.0009 


1.0 




7234l/7677 


0.033 


± 


0.002 [Mo71] 


0.0290 


± 


0.0006 


[Wi80] 


0.0293 ± 0.0011 


1.9 




73019/7677 


n nnn i *? 
U.UUUlo 


_|_ 
It 


U.UUUUo [VVloUJ 










n nnni -U n nnni 
U.UUUl it U.UUU1 






7total/7677 
















0.0334 ± 0.0014 




34 Ar 34 S 


746 1/7666 


0.28 


± 


0.16 [Mo71] 


0.365 


± 


0.036 


[Ha74a] 


0.361 ± 0.035 


1.0 




72580/7666 


0.38 


± 


0.09 [Mo71] 


0.345 


± 


0.01 


[Ha74a] 


0.345 ± 0.010 


1.0 




73129/7666 


0.67 


± 


0.08 [Mo71] 


0.521 


± 


0.012 


[Ha74a] 


0.524 ± 0.022 


1.8 




7total/7666 
















1.231 ± 0.043 




42 Ti 42g c 


72223/7611 


0.012 


± 


0.004 [Ga69] 










0.012 ± 0.004 






7total/7611 


0.023 


± 


0.012 [Ga69,En90] 










0.023 ± 0.012 





a 7-ray intensities are denoted by fE, where E is the 7-ray energy in keV. 



TABLE VI: References for which the original decay-energy results have been updated to incorporate the most recent calibration 
standards. (See Table IVTHl for the correlation between the alphabetical reference code used in this table and the actual reference 
numbers.) 



References (parent nucleus)* 7 



Update procedure 



• Bo64( ls Nc), Ba84( 10 C), Br94( 26 Al r ") 
Ba98( 10 C), Ha98( 3s K m ), To03 ( 14 0) 

• Ry73a( 34 Cl), Sq76( 46 V), Ba77c ( 34 C1) 
Wh77( 14 0) 

• Pr67( 18 Nc) 
. Ja78 ( 38 K m ) 

• Bu79( 3s K m ) 

« Bu61 ( 14 0), Ba62 ( 14 0) 

• Ha74d ( 34 C1) 

• Ba88 and Ba89 ( 10 C) 

• Ki89 ( 42 Sc) 

. Ha74c( 22 Mg, 26 Si, 30 S, 34 Ar) 



• Wc have converted all original (p.n) threshold measurements to Q-values 
using the most recent mass excesses [Au03]. 

• These (p,n) threshold measurements have been adjusted to reflect recent 
calibration a-encrgics [Ry91] before being converted to Q-valucs. 

• Before conversion to a Q-value, this (p,n) threshold was adjusted to reflect a 
new value for the 7 Li(p,n) threshold [Wh85], which was used as calibration. 

• This (p,n) threshold was measured relative to those for 10 C and 14 0; we have 
adjusted it based on average Q-valucs obtained for those decays in this work. 

• Before conversion to a Q-value, this (p,n) threshold was adjusted to reflect the 
modern value for the 35 Cl(p,n) threshold [Au03], which was used as calibration. 

• These 12 C( 3 Hc,n) threshold measurements have been adjusted for updated 
calibration reactions based on current mass excesses [Au03]. 

• These ( 3 He,t) reaction Q-values were calibrated by the 27 Al( 3 He,t) reaction 
to excited states in 27 Si; they have been revised according to modern mass 
excesses [Au03] and excited-state energies [En98]. 

• These measurements of excitation energies in 10 B have been updated to 
modern 7-ray standards [Hc00]. 

• This 41 Ca(p,7) reaction Q-value was measured relative to that for 40 Ca(p,7); 
wc have slightly revised the result based on modern mass excesses [Au03]. 

• These (p,t) reaction Q-valucs have been adjusted to reflect the current Q- 
value for the ie O(p,t) reaction [Au03], against which they were calibrated. 



a These references all appear in Table J] under the appropriate parent nucleus. 
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TABLE VII: References from which some or all results have been rejected even though their quoted uncertainties qualified 
them for inclusion. (See Table IVlHl for the correlation between the alphabetical reference code used in this table and the actual 
reference numbers.) 



References (parent nucleus) 



Reason for rejection 



1. Dccay-cncrgics: 

• Pa72 ( 30 S) 

« No74 ( 22 Mg) 
. Ro74 ( 10 C) 

* Ba77b ( 10 C) 



• Wh81 and Ba98 ( 14 0) 

2. Half-lives: 

. Ja60 ( 26 Al m ), He61 ( 14 0), Ba62 ( 14 0), 
Pr63( 14 0), Fr65b( 42 Sc, 46 V, 50 Mn) 
Si72( 14 0) 

. Ha72a( 26 Al m , 34 C1, 38 K"\ 42 Sc) 
. Ro74 ( 10 C) 

• Ch84 ( 38 K m ) 

3. Branching-ratios: 

• Fr63 ( 26 Si) 



• No calibration is given for the measured (p.t) reaction Q-values; update 
is clearly required but none is possible. 

• Calibration reaction Q-values have changed but calibration process is too 
complex to update. 

• P.H. Barker (co-author) later considered that inadequate attention had 
been paid to target surface purity [Ba84]. 

• P.H. Barker (co-author) later stated [Ba84] that the (p,t) reaction Q-value 
could not be updated to incorporate modern calibration standards. 

• Most of the results in this reference disagree significantly with more recent 
and accurate measurements. Our justification for rejection is presented in 
more detail in the text. 

• The result in [Wh81] was updated in [Ba98] but then eventually withdrawn 
by P.H. Barker (co-author) in [To03]. 

• Quoted uncertainties are too small, and results likely biased, in light of 
statistical difficulties more recently understood (see [Fr69a]). In particular, 
"maximum-likelihood" analysis was not used. 

• All four quoted half-lives are systematically higher than more recent and 
accurate measurements. 

• P.H. Barker (co-author) later considered that pile-up had been 
inadequately accounted for [Ba90]. 

• "Maximum-likelihood 1 ' analysis was not used. 

• Numerous impurities present; result is obviously wrong. 



sets of measurements are simply related to one another by the excitation energy of the isomeric state in the parent. 
In Tabic [I] the set of measurements for the ground-state Qsc-vahie and for the excitation energy of the isomeric state 
appear in separate rows, each with its identifying property given in column 3 and its weighted average appearing 
in column 7. In the row below, the average value given in column 7 for the superallowed transition is the weighted 
average not only of the direct superallowed Qec-vahie measurements in that row, but also of the result derived from 
the two preceeding rows. Note that in all cases the Qsc-vahie for the superallowed transition appears in bold-face 
type. 

For some decays that lead to radioactive daughter nuclei, there is no direct measurement of the Q^c-value for the 
superallowed transition or the one that exists is rather imprecise. In these cases the Qec-value must depend on the 
measured mass excesses of the parent and daughter nuclei, together with the excitation energy of the analog + state 
in the daughter. Each of these properties is identified in column 3 of Table U with the individual measurements of 
that property, their weighted average and a scale factor appearing in columns to the right. The average Qsc-vahie 
listed for the corresponding superallowed transition is obtained from these separate averages. If a direct measurement 
of the superallowed Q^c-value exists, then it is also included in the final average. 

As in our previous survey we have not used the 2003 Mass tables to derive the Q_EC-values of interest. 
Our approach is to include all pertinent measurements for each property; typically, only a subset of the available 
data is included as input to the mass tables. Furthermore, we have examined each reference in detail and either 
accepted the result, updated it to modern calibration standards or rejected it for cause. The updating procedures are 
outlined, reference by reference, in Table IVTl and the rejected results are similarly documented in Table IVTT1 With a 
comparatively small data set, we could afford to pay the kind of individual attention that is impossible when one is 
considering all nuclear masses. 

The half- life data appear in Tabic ITTT1 in similar format to Table [IJ For obvious reasons, half-life measurements do 
not lend themselves to being updated. Consequently, a number of mostly pre-1970 measurements have been rejected 
because they were not analyzed with the "maximum-likelihood" method. The importance of using this technique for 
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TABLE VIII: Reference key, relating alphabetical reference codes used in Tables ITl VIII to the actual reference numbers. 



Table 


Reference 


Table 


Reference 


Table 


Reference 


Table 


Reference 


Table 


Reference 


Table 


Reference 


code 


number 


code 


number 


code 


number 


code 


number 


code 


number 


code 


number 


Ad83 


B 


Aj88 


[10] 


Aj91 


[ii] 


A169 


[12] 


A172 


[13] 


A175 


[14] 


A177 


[15] 


A178 


[16] 


A182 


[17] 


An70 


[18] 


Au03 


[19] 


Az74 


[20] 


Az75 


[21] 


Ba62 


[22] 


Ba67 


[23J 


Ba77a 


[24] 


Ba77b 


[25] 


Ba77c 


[26] 


Ba84 


[27] 


Ba88 


[28] 


Ba89 


[29J 


Ba90 


[30] 


Ba98 


[31] 


BaOO 


[32] 


BaOl 


[33] 


Ba04 


[34] 


Ba06 


[35] 


Ba08 


[36] 


Be68 


[37] 


Be78 


[38] 


Be85 


[39] 


Be08 


[40] 


Bi03 


[41] 


B104a 


[42] 


B104b 


[43] 


Bo64 


[44] 


Br94 


[45] 


Bu61 


[46] 


Bu79 


[47J 


Bu88 


[48] 


Bu06 


[49] 


Ca05 


[50] 


Ch84 


[51] 


C173 


[52] 


Da80 


[53] 


Da85 


[54] 


De69 


[55] 


De78 


[56] 


Dr75 


[57] 


En90 


[58] 


En98 


[59] 


Er06a 


[60] 


Er06b 


[61] 


Er08 


[62] 


Fa09 


[63] 


Fi08 


[64] 


Fr63 


[65] 


Fr65b 


[66] 


Fr69a 


[67] 


Fr75 


[68] 


Fu99 


[(59] 


Ga69 


[70] 


GaOl 


[71] 


Ge07 


[72] 


Ge08 


[73] 


Gi72 


[74] 


Go72 


[75] 


Gr07 


[76] 


Gr08 


[77] 


Ha67 


[78] 


Ha68 


[79] 


Ha72a 


[80] 


Ha74a 


[81] 


Ha74b 


[82] 


Ha74c 


[83] 


Ha74d 


[84] 


Ha75 


[85] 


Ha94 


[86] 


Ha98 


[87] 


Ha02 


[88] 


Ha03 


[89] 


He61 


[90] 


He81 


[91] 


He82 


[92] 


HcOO 


[93] 


HeOl 


[94] 


He02 


[95] 


Ho64 


[96] 


Ho74 


[97] 


Hu82 


[98] 


Hy03 


[99] 


Hy05 


[100] 


Ia06 


[101] 


Ia08 


[102] 


In77 


[103] 


Is80 


[104] 


Ja60 


[105] 


Ja78 


[106] 


Je07 


[107] 


Ji02 


[108] 


Ka68 


[109] 


Ka69 


[110] 


Ke07 


[m] 


Ki89 


[112] 


Ki91 


[113] 


Ko83 


[114] 


Ko87 


[115] 


Ko97a 


[116] 


Ko97b 


[117] 


Kr91 


[118] 


Le08 


[119] 


Li94 


[120] 


Ma94 


[121] 


Ma08 


[122] 


Mc67 


[123] 


Mi67 


[124] 


Mo71 


[125] 


Mu04 


[126] 


Na91 


[127] 


Ni69 


[1281 


No74 


[129] 


OiOl 


[130] 


Pa72 


[131] 


Pa05 


[132] 


Pi03 


[133] 


Pr67 


[134] 


Pr90 


[135] 


Ra83 


[136] 


Re85 


[137] 


Ri07 


[138] 


Ro70 


[139] 


Ro72 


[140] 


Ro74 


[141] 


Ro75 


[142] 


Ro06 


[143] 


Ry73a 


[144] 


Ry91 


[145] 


Sa80 


[146] 


Sa95 


[147] 


Sa04 


[148] 


Sa05 


[149] 


Sc05 


[150] 


Sc07 


[151] 


Se73 


[152] 


Se05 


[153] 


Sh55 


[154] 


Si66 


[155] 


Si72 


[156] 


Sq75 


[157] 


Sq76 


[158] 


Ti95 


[159] 


To03 


[160] 


To05 


[161] 


Vo77 


[162] 


Wa83 


[163] 


Wa92 


[164] 


We68 


[165] 


Wh77 


[166] 


Wh81 


[167] 


Wh85 


[168] 


Wi76 


[169] 


Wi78 


[170] 


Wi80 


[171] 


Zi72 


[172J 


Zi87 


[173] 















precision measurements was not recognized until that time [67[ and, without access to the primary data, there is no 
way a new analysis can be applied retroactively. All rejected half-life measurements arc also documented in Table IVlTl 

Finally, the branching-ratio measurements are presented in Table [TVl The decays of the T z = parents are the most 
straightforward since, in all these cases, the superallowed branch accounts for >99.5% of the total decay strength. 
Thus, even imprecise measurements of the weak non-superallowed branches can be subtracted from 100% to yield the 
superallowed branching ratio with good precision. For the hi gher- Z parents of this type, particularly 62 Ga and heavier, 
it has been shown theoretically [88j and experimentally [64 , 1 1331 ] that numerous very- weak Gamow- Teller transitions 
occur, which, in total, can carry significant decay strength. Where such unobserved transitions are expected to exist 
and have not already been accounted for in the quoted references, we have used a combination of experiment and 
theory to account for the unobserved strength, with uncertainties being adjusted accordingly. 

The branching ratios for decays from T z — — 1 parents are much more challenging to determine, since the super- 
allowed branch is usually one of several strong branches - with the notable exception of 14 - and, in two of the 
measured cases, it actually has a branching ratio of less than 10%. The decays of 18 Ne, 26 Si, 30 S, 34 Ar and 42 Ti thus 
required special treatment. In each case, the absolute branching ratio for a single /3-transition has been measured. 
The branching ratios for other /3-transitions then had to be determined from the relative intensities of /3-delayed 7 
rays in the daughter. The relevant 7-ray intensity measurements appear in Table [Vl with their averages then being 
used to determine the superallowed branching-ratio averages shown in bold type in Tabic IIVI These cases are also 
flagged with a table footnote. 

III. THE Tt VALUES 

With the input data now settled, we can proceed to derive the ft values for the 20 superallowed transitions included 
in the tables. In our last survey @, we described and used a new computer code for calculating the statistical rate 
function /, which surpassed the precision then being obtained with measurements of Qec- Since then, with the advent 
of Penning-trap mass measurements, experimental uncertainties have shrunk even further. The level of precision 
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TABLE IX: Derived results for superallowed Fermi beta decays. 



Parent 




1 t, G 


Partial half-life 










nucleus 


J 


<%) 


t(ms) 


ft(s) 


S'„ (%) 


S„ - Satq (%) 




T z = -1: 
io c 


2.3004 ± 0.0012 


0.297 


1322300 ± 1800 


3041.7 ± 4.3 


1.679 ± 0.004 


0.520 ± 0.039 


3076.7 ± 4.6 


14 Q 


42.772 ± 0.023 


0.088 


71127 ± 51 


3042.3 ± 2.7 


1.543 ± 0.008 


0.575 ± 0.056 


3071.5 ± 3.3 


18 No 


134.47 ± 0.15 


0.081 


21660 ± 590 


2912 ± 79 


1.506 ± 0.012 


0.855 ± 0.052 


2931 ± 80 


22 Mg 


418.39 ± 0.17 


0.069 


7295 ± 17 


3052.0 ± 7.2 


1.466 ± 0.017 


0.605 ± 0.030 


3078.0 ± 7.4 


26 Si 


1029.4 ± 2.2 


0.064 


2954 ± 23 


3041 ± 24 


1.438 ± 0.023 


0.650 ± 0.034 


3065 ± 25 




1966.9 ± 8.0 


0.066 


1524 ± 21 


2998 ± 44 


1.423 ± 0.029 


1.040 ± 0.032 


3009 ± 44 


34 Ar 


3414.5 ± 1.5 


0.069 


894.0 ± 2.4 


3052.7 ± 8.2 


1.412 ± 0.035 


0.845 ± 0.058 


3069.6 ± 8.5 


38 Ca 


5327.2 ± 1.8 


0.075 






1.414 ± 0.042 


0.940 ± 0.072 




42 Ti 


7040 ± 30 


0.088 


470 ± 160 


3300 ± 1100 


1.428 ± 0.050 


1.170 ± 0.080 


3300 ± 1100 



T z = 0: 



26 M m 


478.237 ± 0.038 


0.082 


6350.2 ± 1.9 


3036.9 ± 0.9 


1.478 ± 0.020 


0.305 ± 0.027 


3072.4 ± 1.4 


34 Cl 


1995.96 ± 0.47 


0.080 


1527.77^;« 


3049. 4+ \ I 


1.443 ± 0.032 


0.735 ± 0.048 


3070.6 ± 2.1 


38 K m 


3297.88 ± 0.34 


0.085 


925.42 ± 0.28 


3051.9 ± 1.0 


1.440 ± 0.039 


0.755 ± 0.060 


3072.5 ± 2.4 


42 Sc 


4472.24 ± 1.15 


0.099 


681.43 ± 0.26 


3047.6 ± 1.4 


1.453 ± 0.047 


0.630 ± 0.059 


3072.4 ± 2.7 




7209.47 ± 0.90 


0.101 


422.99 ± 0.11 


3049.5 ± 0.9 


1.445 ± 0.054 


0.655 ± 0.063 


3073.3 ± 2.7 


50 Mn 


10745.97 ± 0.51 


0.107 


283.68 ± 0.11 


3048.4 ± 1.2 


1.444 ± 0.062 


0.695 ± 0.055 


3070.9 ± 2.8 


64 Co 


15766.6 ± 2.9 


0.111 


193 4Q<;+ ' 063 


3050. 


1.443 ± 0.071 


0.805 ± 0.068 


3069.9+3 2 


62 Ga 


26400.2 ± 8.3 


0.137 


116.442 ± 0.042 


3074.1 ±1.5 


1.459 ± 0.087 


1.52 ± 0.21 


3071.5 ± 7.2 


66 As 


32125 ± 470 


0.155 






1.468 ± 0.095 


1.62 ± 0.40 




70 Br 


38600 ± 3600 


0.175 






1.49 ± 0.11 


1.69 ± 0.25 




74 Rb 


47300 ± 110 


0.194 


65.227 ±0.078 


3084.9 ± 7.8 


1.50 ± 0.12 


1.71 ± 0.31 


3078 ± 13 



Average (best 13), Ft 3072.08 ± 0.79 

X 2 /v 0.28 



possible has currently reached ^0.001%, at least for the Qec values of 50 Mn and 54 Co, so it is now necessary to 
include in the / calculation a provision for atomic excitation of the daughter nucleus if the calculation is to continue 
to match the precision of the input data. Our method for including this effect is described in Appendix [XJ and we 
also present there, in Table IXl a comparison of / values both with and without this small correction. It can be seen 
that the effect of the correction is comparable to a shift of 0.001-0.004% in the Qec value, an amount significant 
enough to warrant its inclusion in future. Our final / values are recorded in the second column of Table IIXI They 
were evaluated using our updated code and the Qec values with their uncertainties from column 7 of Table [IJ 

The third column of Table HXl lists (as percentages) the electron-capture fraction, Pec, calculated for each of the 
20 superallowed transitions. The method of calculation was described in our last survey [f|, to which the reader is 
referred for more details. The partial half- life, t, for each transition is then obtained from its total half-life, ti/ 2 , 
branching ratio, R, and electron-capture fraction according to the following formula: 

t= tj g(l + PEc). (2) 

The resultant values for the partial half-lives and the corresponding ft values appear in columns 4 and 5 of the table. 
To obtain the Ft from each ft value, we use Eq. ^ to apply the small transition-dependent correction terms, 
8ns and 8c- We take the values of these terms from our recent re-evaluation of the corrections to superallowed 
beta decay 0]. The first term, S' R , which is listed in column 6 of Table HXl is taken from Tabic V in Ref. J3- The 
two nuclear structure-dependent terms, combined in the form (Sc — Sns), are listed in column 7. In Ref. [7|, 8c is 
expressed as the sum of Sci and 5c2, the former being listed in Table III of that reference and the latter in Table 
II; Sns is taken from Table VI of the same reference. Finally, the resulting Tt values are listed in column 8 of our 
Table El 

Both the uncorrected ft values and the fully corrected Tt values are plotted in Fig. [3] for the 13 most precisely 
measured transitions. The differences between the former, in the top panel, and the latter, in the bottom panel, 
illustrate the effects of our including the correction terms. It is also worth remarking that the values of 5' R are very 
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FIG. 3: In the top panel are plotted the uncorrected experimental ft values as a function of the charge on the daughter nucleus. 
In the bottom panel, the corresponding Tt values are given; they differ from the ft values by the inclusion of the correction 
terms S' R , 8ns and 8c- The horizontal grey band in the bottom panel gives one standard deviation around the average Tt 
value. 



nearly the same for 11 of the 13 cases plotted: only 10 C and 1 O have slightly higher values. Thus, most of the 
differences between the two panels of the figure are due to the effects of the nuclear structure-dependent terms, 5ns 
and 8c- 



A. CVC test 



There are now 13 superallowed transitions whose Tt values have uncertainties less than ±0.4%, with the best case, 
26 Al m , being known an order of magnitude better than that. These data are sufficient to provide a very demanding 
test of the CVC assertion that the Tt values should be constant for all nuclear superallowed transitions of this type. 
The data in column 8 of Tabic HXl clearly satisfy the test, the weighted average of the 13 most precise results being 

Ti= 3072.08 ± 0.79s, (3) 

with a corresponding chi-square per degree of freedom of \ 2 l v — 0.28. That these 13 Tt values form a consistent set 
is also clearly evident from the bottom panel of Fig. [3] Since Tt is proportional to the square of the vector coupling 
constant, G v , then Eq. can be said to confirm the constancy of G v - and to verify this key component of the CVC 
hypothesis - at the level of 1.3 x 10~ 4 . 

Compared with the results of our last survey § , the value of Tt in Eq. ([3]) is somewhat lower but carries a similar 
uncertainty. However, the new analysis is more demanding since, for the first time, it includes the 62 Ga transition, 
which has been very significantly improved in the last four years. This effectively increases the span of masses over 
which the CVC test is being applied; yet even with this addition, the \ 2 l v is actually lower than it was previously. 
The small reduction in the central value of Tt is within the uncertainty of the previous value; it has arisen principally 
from recent changes in the nuclear structure-dependent correction terms [7] . 



13 




,4,3 2 6 A|m =4 C| 3» K m 4 2 g c My 

Parent nucleus 



Mn "Co 



FIG. 4: Summary histogram of the fractional uncertainties attributable to each experimental and theoretical input factor that 
contributes to the final Tt values for the "traditional nine" superallowed transitions. 



B. Tt value error budgets 

We show the contributing factors to the individual Tt- value uncertainties in Fig. [J] for the "traditional nine" cases 
and in Fig. [5] for the remaining eleven. For most of the cases that contribute to the CVC test - 26 Al m to 54 Co in 
Fig. 2] as well as 62 Ga and 74 Rb in Fig. [5]- the theoretical uncertainties are greater than, or comparable to, the 
experimental ones. In these cases, the nuclear-structure-dependent correction, 5c — 8ns, contributes an uncertainty 
of 3-7 parts in 10 4 to all Tt values between 26 Al m and 54 Co but jumps up to 20-30 parts in 10 4 for 62 Ga and 74 Rb 
because of nuclear-model ambiguities. For its part, the nucleus-dependent radiative correction, S' Rl has an uncertainty 
that starts very small but grows smoothly with Z 2 . This is because the contribution to 5' R from order Z 2 a 3 has only 
been estimated from its leading logarithm jl76| and the magnitude of this estimate has been taken as the uncertainty 




FIG. 5: Summary histogram of the fractional uncertainties attributable to each experimental and theoretical input factor that 
contributes to the final Tt values for the eleven other superallowed transitions. Where the error is shown as exceeding 60 parts 
in 10 4 , no useful experimental measurement has been made. 
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in 5' R . As a result, though, for 50 Mn and 54 Co it becomes the leading uncertainty, indicating that a closer look at the 
order Z 2 a 3 contribution to S' R would certainly now be worthwhile. 

For all the transitions from T z =0 parent nuclei, the experimental branching ratios are > 99% and have very small 
associated uncertainties with the exception of 54 Co, which has a 3 x 10 -4 fractional uncertainty attributed to its 
branching ratio, and 74 Rb, which has 10 x 10 -4 . In both cases, this is because they are predicted to have weak 
Gamow- Teller branches that have not yet been observed. We have used an estimate of the strength of the missing 
branches, taken from a shell-model calculation [88| . to assign an uncertainty to the superallowed branching ratio. 
Numerous weak Gamow- Teller branches become an increasingly significant issue for the heavier-mass nuclei with 
A > 62, where they present a major experimental challenge if they are to be fully characterized. Only in the case of 
62 Ga has this been accomplished so far. 

For the decays of 10 C and 14 0, and for all the decays depicted in Fig.[S]except for 62 Ga and 74 Rb, the predominant 
uncertainties are experimental in origin. Many of the experimental branching ratios, and some of the Q-values and 
half-lives have yet to be measured precisely for the cases in Fig. [5l but recent advances in experimental techniques 
have been improving this situation and are likely to improve it even more within the next few years. 



C. Accounting for systematic uncertainties 

So far, we have dealt with the inter-nuclear behavior of Ft-v&hies, examining their constancy as a test of CVC. 
With that test passed at high precision, we are now in a position to use the average Tt-value obtained from these 
concordant nuclear data to go beyond nuclei, obtaining first the vector coupling constant (see Eq. {T])) and then the 
V u d matrix element. Before doing so, however, we must address one more possible source of uncertainty. Though the 
average J-t value given in Eq. ([3]) includes a full assessment of the uncertainties attributable to experiment and to 
the particular calculations used to obtain the correction terms, it does not incorporate any provision for a common 
systematic error that could arise from the type of calculation chosen to model the nuclear-structure effects. To discuss 
this, we divide the problem into two parts: the accuracy of the model as an approximation to the formally complete 
treatment; and the possible existance of systematic uncertainties within the model. 



1. The model approximation 

Very recently Miller and Schwenk [l77| have explored the formally complete approach to isospin-symmetry breaking. 
Their starting point is to define the Fermi matrix element as 

M F = (f\r+\i) = J2(f\°ib a \i) = ^</|4kX^aK>, (4) 

where aj, creates a neutron and b a annihilates a proton in state a. Here i) and |/) are the exact state vectors for the 
full Hamiltonian. If this Hamiltonian commutes with the isospin operators, then \i) and |/} are exact isospin analogs 
of each other, (7r|6 a |z) = (/laJjTr)* and the symmetry-limit matrix element is 

Mo=£|(/kk)| 2 . (5) 

Ot,7T 

If isospin is not an exact symmetry, then \i) and |/) are not isospin analogs and a correction to Mq needs to be 
evaluated. This is the isospin-symmetry-breaking correction, 5c , we seek to determine. It is defined by 

Ml = Ml (1 - 6c) • (6) 

Ideally, to obtain 5c one would compute Eq. using the shell model, and introduce Coulomb and other charge- 
dependent terms into the shell-model Hamiltonian. However, because the Coulomb force is long range, the shell-model 
space would have to be huge to include all the potential states that the Coulomb interaction might connect with. 
Currently this is not a practical proposition. 

To proceed with a manageable calculation, we have developed a model approach [^ Jl78l[i~79T | in which 5c is divided 
into two parts: 

6c = 6ci + 5c2- (7) 

For 5ci, we compute 

^(7|ai| 7 r)( 7 r|6 Q |7)=Af (l- ( 5 cl ) 1/2 , (8) 

a,7r 
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where \t) and |/} are not the exact eigenstates that appear in Eq. ([4]), but are the shell- model eigenstates of an effective 
Hamiltonian (including charge-dependent terms) evaluated in a modest-sized shell-model space. Since this space does 
not allow for nodal mixing, we correct for that limitation by computing the second component, 5c2, obtained from 

^KMWpr^Moa-fe) 172 , (9) 

CK,7T 

where each is a radial overlap integral of proton and neutron radial functions. We justify the efficacy of this second 
term by the following arguments: If the radial functions were identical, then 5c2 would vanish as it should. Otherwise 
the proton radial functions, u p (r), could be expanded in terms of a complete set of neutron functions, u^(r), including 
all possible radial nodes, N: 

v?{r) =J2<*NU n N (r)- (10) 

N 

The isospin-symmctry breaking correction, 5c2, could then be expressed in terms of a^, which from perturbation 
theory could itself be written in terms of matrix elements of the Coulomb interaction. This would be equivalent to 
the nodal mixing included in Eq. ([4]) but left out of the calculation of Sci in Eq. ([8]). The idea is that Sci is the result 
of a tractable shell-model calculation that does not include any nodal mixing, while 6c2 then corrects for the nodal 
mixing that would be present if the shell-model spac e were larg er. 

Clearly our charge-dependent correction terms 0, Il78l 1179} are based on a model, and required approximations 
to make the computation possible. Since no one has yet made a complete calculation without approximations, it is 
impossible to be definitive about any systematic errors that might be introduced by our methods. Only for the lig htest 
superallowed emitter, 10 C, has it been possible so far even to come close to an exact treatment. Caurier et al. |l80l | 
have reported a large no-core shell- model calculation for that system but, even though they were able to extend their 
basis states up to 8hco, their calculated 5c still had not converged to a stable value. However they used their results 
together with perturbation theory to estimate that the full value of 5c should be about 0.19%. This result, which in 
effect used Eq. (fj| and did not split Sc into two parts, agrees completely with our calculated value for 5c — <5ci + ^C2 
of 0.18(2)% (see Table VII in Ref. |7|). This agreement certainly supports the validity of our model. 

Furthermore, it must be noted that our model approach has allowed us to use well-established shell-model and related 
parameters, which were determined from experimental data that are completely independent of the superallowed ft 
values. As is clearly evident from Fig. [31 these calculated corrections do a remarkablejob in converting widely scattered 
ft values into a consistent set of Tt values. Not only that but, as shown in Ref. Q, they also closely reproduce the 
measured results for isospin-forbidden + ^ + (3 transitions in all nuclei for which the shell- model calculation is well 
specified. (This is not the case for 62 Ga.) Of course, although these successes demonstrate that our calculated 5c 
values correctly reproduce the nucleus-to- nucleus variations observed by experiment, they cannot rule out a constant 
shift in the corrections for all nuclei. Even so, it seems highly unlikely that a faulty approximation could lead to 
relative results that are correct in every detail, while being consistently wrong - and by the same constant amount - 
in the absolute values for each and every case. 

Under the circumstances, we see no justification at this time to assign any additional systematic error to account 
for possible inadequacies of the model we use to calculate the charge-dependent correction terms. 



2. Systematic uncertainty within the model 

As introduced in Eq. ([7]) the isospin-symmetry breaking correction, 5c, is separated into two pieces: Sci comes 
from configuration mixing in a modest-sized shell- model calculation with charge-dependent interactions, while 5c2 
involves radial overlap integrals. The calculation of 5ci, the smaller of the two terms, requires a reliable shell- model 
description of the nuclei involved but it can also be further constrained by independent experimental data: for example, 
the measured coefficients of the corresponding isobaric multiplet mass equation (IMME). We take our corrections from 
Ref. [7J where the uncertainties attached to the calculated values of 5ci for the twenty cases of interest here already 
include ample provision for differences between competing shell-model parameterizations as well as for experimental 
uncertainties on the IMME coefficients used as constraints. 

The values of the radial-overlap term, 5c2, which we use as input to Table lIXl were also taken from our recently 
published calculations Q- Those calculations used radial wave functions derived from a Saxon- Woods potential 
with either the well depth or one of the surface terms in the potential adjusted so that the binding energy of each 
computed eigenfunction matched the corresponding measured separation energy. The quoted uncertainties included 
provision for any variations in the results depending on which parameter was used in the adjustment. However no 
provision was included for possible differences that might occur if another method entirely were used to derive the 
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radial wave functions. In the past [5J1, we h ave accounted for this uncertainty by comparing our results for 5c2 with 
those of Ormand and Brown jl8ll Il82l Il83l | , who used Hartree-Fock eigenf unctions and obtained consistently smaller 
corrections than those we found with a Saxon- Woods potential. We treated this as a valid source of systematic 
uncertainty and incorporated it by deriving two average Tt values, one for each set of 5c2 calculations, then taking 
the average of the two and assigning a systematic uncertainty equal to half the spread between them. 

This specific comparison is no longer tenable. The Ormand and Brown calculations are in some cases more than two 
decades old: they use smaller shell-model spaces than are now known to be necessary [7| and they are not available 
at all for some of the transitions that we now need to include. To remedy these deficiencies we have undertaken 
our own Hartree-Fock calculations. They are described in detail in Appendix [Bj where Table IB 21 lists the values 
of Sc2 we compute from Hartree-Fock-derived wave functions and compares them with our earlier results from the 
Saxon- Woods potential 0] , the same results that we used to evaluate Tt in Table [iXj Both methods used exactly the 
same shell-model calculations to determine the full parentage of the states involved. 

With these new Hartree-Fock calculations we can now follow a similar procedure to the one we employed with 
the old calculations in our previous survey p|. We begin by substituting the Hartree-Fock Sc2 values for the Saxon- 
Woods ones in deriving the 5c values used in Table Hxl When we do this the Tt-value result becomes 3071.60±0.89 
with x 2 /V=0.98. This normalized chi-square is three times worse than the one we obtained in Table ITXl with the 
Saxon- Woods corrections, which arguably could justify our rejecting the Hartree-Fock results outright. However, to 
be safe, we proceed as before and take the average of the Hartree-Fock and Saxon- Woods results, adding a systematic 
uncertainty equal to half the spread between the two results. Thus, we obtain 

Tt = 3071.87 ± 0.79 sta t ± 0.27 syst s 

= 3071.87 ±0.83 s, (11) 

where on the second line the two uncertainties have been added in quadrature. Our new systematic adjustment 
amounts to only 0.27s, much smaller and of opposite sign to the 0.90s correction applied previously [f|. 

It is the value for Tt in Eq. fTTj) that we carry forward to subsequent sections where we obtain V u d and test the 
unitarity of the CKM matrix. 



IV. THE IMPACT ON WEAK-INTERACTION PHYSICS 



A. The Value of V ud 



With a mutually consistent set of Tt values, we can now use the adjusted average value in Eq. (JTTJ) to determine 
the vector coupling constant, G v , from Eq. The value of G v itself is of little interest but, together with the weak 
interaction constant for the purely leptonic muon decay, G F , it yields the much more interesting up-down element 
of the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix. The basic relationship is V ud = G V /G F , which in 
terms of Tt becomes 

T/2 = K 

ud 2Gl(l + Al)Tt' V 1 

where is the nucleus-independent radiative correction. This correction has recently been carefully re-examined by 
Marciano and Sirlin Q , who very substantially reduced its uncertainty. Expressing their new result in a way that is 
consistent with the definition of our other correction terms, we obtain (see Eq. (41) in Ref. 0]) 

= (2.361 ± 0.038)%. (13) 

Using the Particle Data Group (PDG) [l74j value for the weak interaction coupling constant from muon decay of 
G F /(hc) 3 = (1.16637 ±0.00001) x 10~ 5 GeV" 2 , we obtain from Eq. H2]) the result 

\V ud \ 2 = 0.94916 ± 0.00044. (14) 

Note that the total uncertainty here is almost entirely due to the uncertainties contributed by the theoretical correc- 
tions. By far the largest contribution, 0.00035, arises from the uncertainty in A^; 0.00020 comes from the nuclear- 
structure-dependent corrections 5c — 5ns and 0.00004 is attributable to 5' R . Only 0.00016 can be considered to be 
experimental in origin. 

The corresponding value of V u d is 

\V ud \ = 0.97425 ± 0.00022, (15) 
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FIG. 6: Values of V u d as determined from superallowed + — > + /3 decays plotted as a function of analysis date , span ning the 
past two decades. In order, from the earliest date to the most recent, the values are taken from Refs. |l84l ]. [l85l | . [f| and 
this work. 

a result which is consistent with, but more precise than, values we have obtained in previous analyses of superallowed 
8 decay. To emphasize the consistency and steady improvement that has characterized the value of V u d as derived 
from nuclear 8 de cay, in Fig . [6] we plot our new result together with V u d values published at various times over the past 
two decades [3, H, Il84l Il85l | Clearly the large number of measurements that contribute to the nuclear determination 
of Vud provides a robust data base, one not subject to sudden shifts as each new measurement appears. 

B. Unitarity of the CKM matrix 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix transforms one set of quark basis states into another: it transforms 
the quark-mass eigenstates into the weak-interaction eigenstates. If both sets are complete and orthonormal, then the 
transformation matrix itself must be unitary. The Standard Model does not prescribe the individual elements of the 
CKM matrix - they must be determined experimentally - but absolutely fundamental to the model is the requirement 
that the matrix be unitary. Whether unitarity is satified in practice can be tested experimentally, the severity of the 
test depending on the precision with which the CKM matrix elements can be determined. 

To date, the most demanding test of CKM unitarity comes from the sum of squares of the top-row elements, 
\V u d\ 2 + \V US \ 2 + IKtf,! 2 , which should equal one. Since \V u d\ 2 constitutes 95% of this sum, the precision on V u d is of 
paramount importance. The value of \V u d\ = 0.97425(22) derived in Sect. IIV Al has a precision of 0.02% which is the 
most precise result so far obtained for this matrix element and is, by more than an order of magnitude, the most 
precisely determined value for any element in the CKM matrix. Alternative methods of obtaining V u d from neutron 
beta decay, V u d = 0.9746(1 8), a nd from pion beta decay, V u d — 0.9749(26) - both values taken from the Particle Data 
Group's 2008 compilation jl74| - are much less precise and have been hampered by experimental difficulties. In the 
case of the neutron, not only is physical containment a problem but the axial-vector contribution to its 8 decay must 
be separated from the vector contribution by a beta-asymmetry measurement; while for pion beta decay a very small 
branching fraction, 0(1O -8 ), must be measured. 

At the time of our last survey [|[ the value of V us was in a state of flux. The 2004 Particle Data Group value 
was \V US \ — 0.2200(26), based mostly on measurements that were at least two decades old, but new results then 
emerging were suggesting a value some two standard deviations higher. In the last f our y ears, th ese n ew resul ts on 
the semi-leptonic decays, Km, of both charged and neutral kaons - from BNL-E865 fl86j . KTeV [l87j], NA48 (l8§ |. 
KLOE [189( | . and IST RA+ [l90| | - have all combined to clarify the situation. Now, current averages by the 2008 
Particle Data Group [174[ and FlaviaNet Q for kaon semi-leptonic branching fractions are based only on recent, 
high-statistics experi ment s, which are also consistent with one another. The best current value, presented at the 
CKM2008 Workshop [l9l| by the FlaviaNet group, is 

/+ (0) | Vus | =0.21673 ±0.00046. (16) 

Here /+(0) is the semi-leptonic decay form factor at zero-momentum transfer. Its value is close to unity. In fact, 
the CVC hypothesis in the exact SU(3) symmetry limit establishes its value to be exactly one, but SU(3) symmetry 
is broken to some extent and a theoretical calculation is required to estimate the departure of /+(0) from unity. 
Curren tly, there are two class es of evaluation: analytic or semi -ana l ytic appr oaches based on chiral perturbation 
theory [Si [HI [HI GH, [l9(| , and those based on lattice QCD [Hll, [HI \l9^ . We will follow the FlaviaNet group 
and adopt the lattice value of /+(0) = 0.9644 ± 0.0049 from the RBC-UKQCD collaboration [113], which yields 




\V US \ = 0.2247 ±0.0012. 



(17) 
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An independent determination of V us can be obtained from the purely leptonic decay of the kaon, the most important 
mode being K + — ► [i^v. If it is considered as a ratio with the leptonic decay of the pion, 7r + — > the hadronic 

uncertainties can be minimized and the result yields the ratio of the CKM matrix elements |Vus|/|Kid|- In the analysis 
of the FlaviaNet group [8J] the current result is 

W f 

' = 0.2760 ±0.0006, (18) 



I Vud | /tt 

where fx and are the kaon and pion decay constants. This ratio of pseudoscalar decay constants has to be obtained 
from theory, for which lattice QCD seems to be the only reli able source. Again, following the FlaviaNet group we 
adopt the lattice result from the MILC-HPQCD collaboration [200j], f K /U = 1.189 ± 0.007, and obtain 



it; 



I ( ..s 



\v ud 



= 0.2321 ±0.0005. (19) 



Thus, we now have three pieces of data - \V u d\ from nuclear decays, Eq. (fTB")) . T4s| from decays, Eq. (fT7| . 
and the ratio 1 14 s | / 1 Vud | from Kn decays, Eq. (fT9|) - from which to determine two parameters, \V u d\ and \V US \- We 
perform a non-linear least squares fit to obtain the result 

\V ud \ = 0.97424(22) \V US \ = 0.22534(93). (20) 

Note that the value of \V ud \ obtained from this fitting procedure has only changed by one unit in the last figure 
compared to Eq. (fT5|) ; and the change in \V US \ compared to Eq. (|17|. though somewhat larger, is still well within the 
quoted uncertainties. 

The third element of the top row of the CKM matrix, V u b, is very small and hardly impacts on the unitarity test 
at all. Its value from the 2008 Particle Data Group compilation jl74| is \V u b\ = (3.93 ± 0.35) x 10~ 3 . Combining this 
number with the ones in Eq. (120(1 we find the sum of the squares of the top-row elements of the CKM matrix to be 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.99995 ± 0.00061, (21) 

a result that shows unitarity to be fully satisfied at the 0.06% level. Only V us and V u d contribute perceptibly to the 
uncertainty and their contributions are almost equal to one another. This may seem surprising since V u d is known to 
much higher precision than V us , but it follows from the fact that |V^d| 2 contributes 95% to the unitarity sum. 



C. Limit on Scalar Interactions 

1. Fundamental scalar current 

In our previous survey || we explained in detail how a scalar current, if it existed, would affect the Tt- value data. 
We demonstrated that its effect on Tt would be approximately proportional to (1/W), the average inverse decay 
energy of each j3 + transition, so its presence would be manifest by Tt values that are not constant as a function of 
Z. Since (1/W) increases monitonically as Z decreases, the largest deviation of Tt from constancy would occur for 
the superallowed transitions from nuclei with the lowest Z, 10 C and 14 0. 

We have now repeated the same analysis on our new survey results. We evaluated the statistical rate functions, 
/, with a shape-correction factor that included the presence of a scalar current via the Fierz interference term, bp, 
which we treated as an adjustable parameter. We then sought the value of &f that minimized x 2 m a least-squares 
fit to the expression Tt — constant. The result we obtained is 

b F = -0.0022 ±0.0026, (22) 

which is consistent with zero, as it was in 2005 In Fig. [7] we illustrate the sensitivity of this analysis by plotting 
the measured Tt values together with the loci of Tt values that would be expected if &f = ±0.004. Obviously, the 
measured Tt values do not exhibit any statistically significant curvature. 

Th e result in Eq. (|2"2")) can also be expressed in terms of the coupling constants that Jackson, Treiman and Wyld 
[20ll | used in writing a general form for the weak-interaction Hamiltonian. Since we are dealing only with Fermi 
superallowed transitions, we can restrict ourselves to scalar and vector couplings, for which that Hamiltonian becomes 
the following: 

H s +v = {ip p ipn)(Cs<f> e <he + C's<i>el$<foe) 
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FIG. 7: Corrected Tt values from Table HXl plotted as a function of the charge on the daughter nucleus, Z. The curved lines 
represent the approximate loci the J-t values would follow if a scalar current existed with bp = ±0.004. 



where we have taken the vector current to be maximally parity violating, as indicated by experiment [202j. The 
complexity of the relationship between bp and the couplings Cs, C' s and Cy depends on what assumptions are made 
about the properties of the scalar current. If we take the most restrictive conditions, that the scalar and vector 
currents are time-reversal invariant (i.e. Cs and Cy are real) and that the scalar current, like the vector current, is 
maximally parity violating (i.e. Cs = C' s ), then we can write 



Cs bp 

= —i- = +0.0011 ± 0.0013. 

Cy 2 



(24) 



This limit from superallowed (3 decay is, by far, the tightest limit on the presence of a scalar current under the 
assumptions stated. 

If we remove the condition that the scalar current be maximally parity violating, then the expression contains two 
unknowns, 



bp 



-2Cy(C s + C s ) 
\Cs\ 2 + \C' s \ 2 
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C's 



(25) 



and cannot be solved individually for Cs/Cy and C' s /Cy. However, the f3-v angular-correlation coefficient, a, for a 
superallowed + ^ + [3 transition provides another independent measure of Cs and Cy. In that case, 
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(26) 



which, together with Eq. ([23]) . can be used set limits on both Cs/Cy and C'g I Cy . Cu rrently, the most precise 
measurement of such a {3-v angular correlation is for the superallowed decay of 38 K |203l |. In this case, what was 
actually measured is a = a/(l + jbpirie/ (W)), where 7 = yl — (aZ) 2 and m e is the mass of th e electron 1 . The 
results in terms of Cs/Cy and C' s /Cy are plotted in Fig. [51 The value of a taken from Ref. [203| leads to the grey 



annulus plotted in the figure, while our result for bp from Eq. ((22)) is responsible for the narrow diagonal band. The 
intersection of these two regions, which is in black, defines the 68% confidence limit (one standard deviation) for 
Cs/Cy and C' s /Cy. It corresponds to the limit 



Cs_ 

Cy 



< 0.065, 



(27) 



and exactly the same limit for C' s /Cy. 



1 Our bp is defined differently from the b used in Ref. |203H . The two are related by bp = 6/7. 
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C s /C v 

FIG. 8: Allowed range of values for Cs/Cv and C' s /Cv as determined from the 38 K f3-v angular-correlation measurement [203| ] 
(grey annulus) and from superallowed + — > + (3 decays (narrow diagonal band). The overlap region is shaded in black. We 
plot the 68% confidence limits. 

The Jackson- Treiman-Wyld Hamiltonian [20l| is a parameterization of one possible extension to the Standard 
Model. The coupling constants Cs and C' s are not prescribed by the authors but are simply parameters that must 
be determined from experiment. One model that actually introduces scalar interactions in a natural way is the 
minimal supersymmetric standard model (MSSM) . A valuable review of low-energy tests of this model has recently 
been published by Ramsey-Musolf and Su jl77l |. In the MSSM a radiative correction to the beta-decay amplitude 
involves box graphs with exchanged supersymmetric sfcrmions. These produce an energy dependence in the beta- 
decay amplitude that shows up in the Fierz interference term. Unfortunately, estimates by Profumo, Ramsey-Musolf 
and Tulin [205j | indicate that the resulting value of &f would be less than 10 -3 , which is an order of magnitude smaller 
than our current experimental upper limit on that quantity. Thus it will likely be some time before superallowed beta 
decay can provide useful constaints for this class of supersymmetric models. 

2. Induced Scalar current 

If we consider only the vec tor p art of the weak interaction, for composite spin-1/2 nucleons the most general form 
of that interaction is written [206l | as 

H v = i> p {gvlv - fu^fiuq v + ifsQii)ip n delvi 1 + l5)<fc e , (28) 

with being the four- momentum transfer, g M = (p p —p n )n- The values of the coupling constants g v (vector), f M 
(weak magnetic) and / s (induced scalar) are prescribed so long as the CVC hypothesis - that the weak vector current 
is just an isospin rotation of the electromagnetic vector current - is correct. In particular, since CVC implies that 
the vector current is divergenceless, it follows that / s should equal zero. An independent argument [207j |. that there 
be no second-class currents in the hadronic weak interaction, also requires / s to vanish. We proved in our previous 
survey [H that the presence of a non-zero / s would manifest itself in exactly the same way as a non-zero Cs- by a 
(1/W) dependence in the Tt- value data. 
In the same manner that we obtained Eqs. (|22[) and (|24p . we determine from our present survey results that 

m e f s /g v = -(0.0011 ± 0.0013). (29) 

This result is a vindication for the CVC hypothesis, which predicts g y = 1 and f s = 0. We confirm this prediction 
at the level of 24 parts in 10 . Our result can also be interpreted as setting a limit on vector second-class currents in 
the semi-leptonic weak interaction. 

D. Limits on extensions to the standard model 



The unitarity sum established in Sect. IIV Bl can be used to set limits on new physics beyond the standard model. 
A list of possible extensions includes, but is certainly not limited to, right-hand currents, extra Z bosons, scalars, 
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supersymmetry, a fourth gen erati on of quarks and exotic muon decay. Marciano surveyed many of these possibilities 



at the CKM2008 Workshop [191j . His general conclusion was that although the CKM unitarity test yields no sign 



of new physics, it does place important constraints on the possibiliti es. In th e cas e of supersymmetric models these 
constraints have been explored by Ramsey-Musolf , Su and Kurylov |204l . 12081 1209| | . In the minimal supersymmetric 
version (MSSM), corrections to low-energy observables arise only via loop effects, while in extensions that allow for 
i?-parity violating (RPV) interactions new tree-level effects appear. In general, the presence of new physics may 
modify low-energy semi-leptonic electroweak observables in two ways: (i) directly, via a new semi-leptonic interaction 
(e.g. right-hand currents, MSSM with RPV interactions), and (ii) indirectly, via loop graphs contributing to the 
radiative correction (e.g. extra Z-bosons, MSSM). In what follows we give one example of each type of modification: 
right-hand currents and extra Z bosons. 



1. Right-hand currents 

In the standard model, parity viola tion i s co nsidered to be maximal. What if this condition were to be relaxed? 
For semi-leptonic transitions, Herczeg [210l l21l| extends the general form of the weak interaction to read 

H se = a LL (V - A)(V — A) + a LR (V - A)(V + A) + a RL (V + A)(V - A) + a RR (V + A)(V + A), (30) 

where, in each term, the first factor represents the lepton currents and the second, the hadron currents. In particular, 
for the vector lepton current, V stands for either fy^^h or <fi e j^(j)3 depending on whether the chirality of the 
neutrino is left-handed, as it is for V — A coupling, or right-handed, as it is for V + A coupling. In the standard 
model, cill = 1, and cil R = a R L = clrr = 0. For Fermi beta decay, only the vector part of the weak hadron current 



y ti 



contributes, so the decay rate is given by the following proportionality 

T/3 oc \a L L + a L R,\ 2 + \a RL + a RR \ 2 

~ \a LL | 2 (l + 2Rea LR + ...), (31) 

where cllr = cllr/o-ll- In the second line of the equation, we have only retained quantities that are first order in the 
(presumably) small quantities cllr, Hrl and Hrr. 

To determine the effect that right-hand currents would have on the value of V ud obtain ed fr om experiment, we 
also need to consider the role of such currents on the purely leptonic muon decay. Herczeg |21C| writes the effective 
Hamiltonian, in analogy to Eq. (|30j) . as 

H e = c LL (V - A)(V -A)+ c LR (V - A)(V + A) + c RL (V + A)(V - A) + c RR (V + A)(V + A). (32) 

The coupling constants in Eqs. (j32|) and (f30|) are related by the CKM matrix elements: 



O-LL — CLLVud a LR — CLR^V^ 

a R L = c RL V^ d a RR = CRRe la V^ d . (33) 

Here V^ d is the element of the CKM matrix for left-handed chirality quarks, and V^ d is for right-handed chirality 
quarks. The phase a is a CP-violating phase in the right-handed CKM matrix. The decay rate for muon decay is 
constructed from an equal mix of vector and axial- vector interactions and is proportional to the following expression 
i: 

oc \c L l\ 2 + \clr\ 2 + \crl\ 2 + \crr\ 2 

= \cll\ 2 {1 + \clr\ 2 + \crl\ 2 + \c R r\ 2 ), (34) 

where Cy = c^/cll- 

If we define \V ud \ 2 expt as being the quantity obtained from the ratio of measured beta- and muon-decay rates, we 
can combine Eqs. (|3"Tj) and ([34| to relate this experimental result to the matrix element |T^| 2 by the relationship 



iy |2 ^ lyL |2 I 1 + a LR\ 2 + \a-RL + QRb| 2 

*ud\expt — t-\ " ud\ -1,1— I? i I— 19 i I— 12 

T M 1 + \clr\ 2 + \c R l\ z + \Crr\ z 

tL i2 



V^\ 2 (l + 2Rea LR ), (35) 



where, in the second line, only corrections to first order in small quantities are retained. If the situation is identical 
for the second (kaon decay) and third (B-meson decay) generations of quarks, with the interaction coupling constants 
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ciij and Cij in H s g and Hi being generation independent, then 

i i 

= l + 2Rea LR . (36) 

In writing the second line we have assumed that the CKM matrix for left-hand chirality quarks is strictly unitary. 
Since the left-hand side of Eq. ([56]) is the experimentality determined unitarity sum, given in Eq. (|2"Tj) of Sect. IIV Bl 
this expression can clearly be used to set a limit on the coupling constant ollr- The result is 

0.99995 ±0.00061 = 1 + 2Rea LR 

Rea LR = -0.00003 ±0.00030, (37) 

which is consistent with no right-hand currents - at least not in the LR sector. 



2. Extra Z bosons 



The existence of neutral gauge bosons, beyond the usual photon and Z boson of the standard SU (2)l x Z7(l) model, 
would impact on the CKM unitarity test. To illustrate this we consider just one of the many models that appear in 
grand unified theories, namely the 50(10) model, whose group breakdown is 

5O(10) -> SU(3) C x SU{2) L x U(l) x U(l) x . (38) 

Here an extra U(l) group is introduced, U(l) x , and its corresponding neutral gauge boson is labelled Z x . The 
existence of such an extra Z boson would impact on the calculation of the electroweak radiative correction. One of 
the important diagrams in the hadron- independent radiative correction, A^, is a WZ -box graph. This graph would 
have to be augmented by an additional WZ x -box graph, whose contribution is of order lax x , where x x — m 2 z /m^, 
the ratio of squared masses of the heavy bosons in the box diagram. If we assume that this correction is common to 
all quark flavours, then the same correction that occurs in the determination of K«;| 2 would also occur for |V^ S | 2 and 
Ktbj 2 . If so, its impact can be in corp orated into the unitarity test. 
Following Marciano and Sirlin [212| . we write 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = (0.99995 ± 0.00061) + A, (39) 

where the numerical value for the experimental unitarity sum is from Eq. I]2ip. and A is a calculated correction due 
to the extra Z boson. If we take the CKM matrix to be exactly unitary in three generations, then Eq. (f3"9"|) can be 
used to set the following one-standard-deviation limits on A: 

- 0.00056 < A < +0.00066. (40) 

Marciano and Sirlin [212j | have computed the contribution of a putative Z x boson to the radiative correction and 
obtained 
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A = 5 x - C, 2 \, (41) 

407rsin 2 <V 3 1 xl x x -V v ' 

where a is the fine-structure constant, Qy/ is the Weinberg angle (sin 2 9w — 0.23), and C x is a coupling constant 
linking the Z x boson to fermions. The normalization has been selected so that C x is unity at the 5O(10) unification 
mass scale. Its value at lower energies has to be estimated, and Marciano and Sirlin use |O x | 2 = \. Noting that the 
correction A is negative, we obtain from the lower limit in Eq. (|40p 

< 0.12. (42) 

x x - 1 

Taking for the W-boson mass, mw — 81 GeV, we arrive at the limit 

m Zx > 430 GeV. (43) 

Impressive though this limit is, somewhat higher limits have been obtained in direct searches at proton and electron 
colliders. The CDF and DO experiments at FermiLab in searches of pp — > e + e~ have placed lower-mass limits (at 95% 
C.L.) on mz x of 822 and 640 GeV respectively, while at CERN the LEP2 experiment on e + e~ — > // (with / signifying 
a fermion) find a lower-mass limit of 673 GeV. These limits are recorded in the survey of Erler and Langacker [2131 ] 
in the 2008 Particle Data Group listings. 
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FIG. 9: Experimental ft values plotted as a function of the charge on the daughter nucleus, Z. Both bands represent the 
quantity Tt/((1 + 5' R )(1 — Sc + 8ns))- The two separate bands distinguish those beta emitters whose parent nuclei have isospin 
T z — — 1 (darker shading) from those with T z — (lighter shading). 

V. CONCLUSIONS 

In our previous survey only four years ago, we remarked on the excellent agreement among the derived Tt values, 
lamented that the results of the unitarity test were still ambiguous, and predicted that the already well-measured ft 
values of the "traditional nine" superallowed decays were unlikely to be improved dramatically in the near future. 
Much has happened since then, not all of it expected. Today, we can say that the excellent Tt-vahie consistency 
remains - or, to be more accurate, it has been restored after Penning-trap Q^c-value measurements, non-existant at 
the time of the last survey, did in fact make important improvements (and changes) in the known ft values, which 
in turn prompted improvements (and changes) in the calculated isospin-symmetry-breaking corrections. At the same 
time, the calculation of the nucleus- independent radiative correction, A^, was improved, leading to a more precise 
result for V u d, and the kaon-decay community mounted a concerted effort, which led to a new and reliable value for 
V us - With these new results, and others, CKM unitarity has now been tested to unprecedented precision . . . and it 
has passed the test with flying colors. 

Furthermore, we have demonstrated in Sec. lIVI how powerful these improved results can be in setting limits on new 
physics beyond the standard model, whether that new physics be a scalar interaction, right-hand currents or extra Z 
bosons. We have seen that tiny uncertainties on the ft values are essential ingredients of a demanding test of CKM 
unitarity, which also leads to tight limits on new physics. The challenge now is: Can those uncertainties be reduced 
still farther? The motivation is as strong as ever: to identify the need for new physics - or to limit the possible 
candidate theories even more definitively. 

We have taken pains throughout this work to pay careful attention to all uncertainties, theoretical and experimental. 
In Sec. IIV Al we detail the various contributions to the uncertainty in |T4d| 2 - Of these, by far the largest is still from 
A^, even though its uncertainty has recently been improved significantly [f|. To improve it more must remain an 
important theoretical goal. 

The next largest contributor to the error budget for \ V u d\ 2 is the nuclear-structure-dependent corrections, (6c — 5ns)- 
Their uncertainties arise both from the input parameters used in their calcul ation - two-body matrix elements in 
the shell-model calculations, experimental uncertainties in charge-radii, etc @, Il79l | - and from possible systematic 
differences between two different methods used for calculating radial wave functions (see Sec. lIII C|) . From a theoretical 
point of view, it would obviously be desirable to have a third completely different calculation, to reinforce the 
assessment of systematic uncertainties. However, in the absence of such a calculation, one must rely on experiment to 
test the accuracy of these calculated corrections. This has become, and should remain, a top priority for experiment. 

The method, which is best described with reference to Fig. [5J is based on the validity of the CVC hypothesis that 
the corrected Tt values for the superallowed + — ► + decays should be constant. In the figure we compare the 
uncorrected measured ft values (points and error bars) with the quantity Tt/{{\ + S' R )(1 — 5c + 6ns)) shown as a 
band, the width of which represents the assigned theory error. The band corresponds to the calculated corrections 
normalized to the data via the measured average Tt value, Tt, taken from Table HXl Thus, although this comparison 
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does not test the absolute values of the correction terms, it does test the collective ability of all three calculated 
correction terms to reproduce the significant variations in ft from one transition to another. In fact, since S' R is 
almost independent of Z when Z > 10, this test really probes directly the effectiveness of the calculated values of 
(S c - S NS )- 

It can be seen that there is remarkable agreement between theory and experiment. In assessing the significance 
of this agreement, it is important to recognize that the calculations of 5c and Sns for Z < 26 are based on well- 
established shell-model wave functions that were further tu ned t o reproduce measured binding energies, charge radii 
and coefficients of the isobaric multiplet mass equation 0, 1 X 791 ) . The origins of the calculated correction terms for 
all cases are completely independent of the superallowed decay data. Thus, the agreement in the figure between 
the measured superallowed data points and the theoretical band is already a powerful validation of the calculated 
corrections used in determining that band. The validation becomes even more convincing when we consider that it 
would require a pathological fault indeed in the theory to allow the observed nucleus- to- nucleus variations in Sc and 
Sns to be reproduced in such detail while failing to obtain the absolute values to comparable precision. As satisfactory 
as the agreement in Fig. [9] is, though, new experiments can still improve the test, making it even more demanding, 
and can ultimately serve to reduce the uncertainty in the nuclcar-structure-dependent corrections even further. 

These new experiments can follow different paths. In the last four years, the biggest impact has come from experi- 
ments that focused on the "traditional nine" superallowed transitions. New Penning-trap Qsc-vahie measurements 
have already been mentioned, but there have been new half-life and branching-ratio measurements as well (see Ta- 
bles [HI] and |IV| . More improvements are still possible as a glance at Fig. Q] reveals. If we accept as a goal that 
experiment should be more than a factor of two more precise than theory, then we see that the Qec values for 
10 C, 14 and 34 C1, the half-lives of 26 AP 1 , 34 C1, 42 Sc and 50 Mn, and the branching ratios for 10 C and 14 can all 
bear improvement. It is also particularly noteworthy that any improvements in the cases of 10 C and 14 will lead 
directly to improvements in the limits on the possible existence of scalar currents. As is evident from Fig. [7] and the 
discussion in Sec. IIV C 1[ it is on these two low-Z superallowed transitions that a scalar current would have the largest 
effect. Unfortunately the branching ratios for both these transitions offer experimental obstacles that have proved 
very difficult to surmount. 

A second experimental path is to expand the number of precisely measured superallowed emitters to include cases 
for which the calculated nuclear-structure-dependent corrections are larger, or show larger variations from nuclide to 
nuclide, than the values applied to the "traditional nine" cases. We argue that if the experimental ft values agree 
with the calculations where the nucleus-to-nucleus variations are large, then that must surely verify the calculations' 
reliability for the nine cases whose corrections are considerably smaller. Already four cases of this type have been 
carefully measured, 22 Mg, 34 Ar, 62 Ga and 74 Rb. They appear to agree well with the calculations although, with 
the exception of 62 Ga, their uncertainties are still five times greater than those for the best known transitions. 
Undoubtedly these uncertainties will be reduced and more cases added in the near future. 

These new cases certainly present serious experimental challenges. The parent nuclei are more exotic than the 
traditional cases, which all have stable daughters, so they are more difficult to produce in pure and statistically 
significant quantities. They also exhibit more complex branching patterns: Each Tz = — 1 parent nucleus decays by 
Gamow- Teller transitions of comparable strength to the superallowed Fermi one, thus requiring the latter's branching 
ratio to be measured directly with high precision. For the Tz = parents with A > 62, each decay includes numerous 
weak Gamow- Teller transitions, which are very difficult to observe individually but which collectively constitute 
nonncgligiblc branching strength. In both regions, these problems are being, or have been overcome, albeit with very 
specialized techniques. The recently published branching-ratio measurement |64ll for 62 Ga is an example of how even 
meticulously detailed spectroscopic studies must be combined with theory [881 ] to ensure that missing transitions are 
properly accounted for in the decays of the heavy Tz — parents. 

There is a further important issue that arises for the superallowed emitters with A > 62: The shell- model calculations 
of the structure-dependent corrections for these nuclei are not solidly based on spectroscopic measurements as they are 
for the lighter nuclei. Such measurements simply do not exist for most N ~ Z nuclei in this mass region. Furthermore, 
charge radii and coefficients for the isobaric multiplet mass equation are not known either and so cannot be used to 
constrain the radial wave functions or "tune" the charge-dependence embedded in the two-body matrix elements. 
As a consequence, the uncertainties assigned to the calculated corrections are very large (see the broad band in this 
mass region in Fig. [9]), considerably reducing the usefulness of these nuclei either in testing the corrections or in 
contributing to the determination of V u d- ft would be very valuable in this context for radioactive-beam facilities to 
direct some attention to determining a wide variety of spectroscopic information in this mass region with a view to 
obtaining a reasonably effective nuclear model, which, among other things, could lead to much improved calculations 
for the correction terms. 

In conclusion, we can assert - as we did four years ago - that world data for superallowed + decays strongly 
support the CVC expectation of an unrenormalized vector coupling constant, and also set a tight limit, consistent 
with zero, on scalar currents. We can now add, though, that CKM unitarity is satisfied to within an uncertainty of 
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0.06%. This reconciliation with unitarity has come about as a result of significant changes in V us ; the value of V u d 
determined from nuclear (3 decay has not varied outside of error bars in twenty years, during which time the size of 
those error bars has been reduced by a factor of five. Finally we have noted that the calculated nuclear-structure 
dependent correction terms have recently been improved and continue to stand up favorably to experimental tests, 
an outcome that must further increase confidence in the nuclear results. 
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APPENDIX A: ATOMIC OVERLAP CORRECTION TO THE STATISTICAL RATE FUNCTION 

The statistical rate function, /, is an integral over phase space, 

pW 

f = J pW(W Q - W) 2 F(Z, W)S(Z, W) dW, (Al) 

where W is the total energy of the electron in electron-rest-mass units; Wq is the maximum value of W; p — (W 2 — l) 1 / 2 
is the momentum of the electron; Z is the atomic number of the daughter nucleus; F(Z, W) is the Fermi function and 
S(Z, W) is the shape-correction function. The details of the calculation of S(Z, W) were given in our previous survey 
Q and will not be repeated here. What we address here is the inclusion for the first time of an additional factor in 
Eq. (|A1[) to account for the mismatch in the initial and final atomic states in the f3 decay. Since the nucleus changes 
charge by one unit in beta decay, the final atomic state does not overlap perfectly with the initial atomic state, an 
effect that leads to a slight inhibition in the beta-decay rate. In the past, this effect has justifiably been considered 
too small to be of practical concern but, with the advent of Penning-trap mass measurements, the experimental 
uncertainties in transition Q-values have been reduced so much that they are now comparable to the effects of the 
imperfect atomic overlap. 
We begin by writing 

f = J pW(W - W) 2 F(Z, W)S{Z, W)r(Z, W) dW, (A2) 

where r(Z, W) is the atomic overlap correction we are seeking. We then follow the method of Bahcall |214j by 
expressing / as a double integral with an energy-conserving delta function: 




where q is the neutrino momentum. We have introduced into this equation an overlap of the initial and final atomic 
electron configurations: \G) is the state vector for the initial neutral atom with (Z + 1) electrons, and \A') is the 
state vector for the final ionized atom with (Z + 1) electrons but only charge Z in the nucleus. There are many such 
possible final states, so a sum over A' is included. 

Of the two energies within the delta function, the first, Ei 1 is the energy of the initial neutral atom in its atomic 
ground state: 

Ei = M z +i{G) = M z+1 + (Z + l)m e - B(G), (A4) 

where Aiz+i{G) is the atomic mass, Mz+i is the nuclear mass, m e is the electron mass (m e = 1 in electron rest-mass 
units) and B(G) is the total electron binding energy in the ground state of the atom. The sign of the latter is chosen 
so that B(G) > 0. The second energy, Ef, is that of the final state, which is composed of an ionized atom still with 
(Z + 1) atomic electrons in an excited configuration plus an emitted beta-decay positron and an emitted neutrino: 



E f = Mzl^A') +W + q = M Z + (Z+ l)m e - B{A') + W + q, 



(A5) 
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where Mz+i( A ') ^ s the atomic mass of a negatively ionized atom (superscript —1 denotes ionization) of (Z + 1) 
electrons in configuration A', Mz is the nuclear mass for the final nucleus in the beta decay, and B(A') is the total 
electron binding energy for the ionized atom. Thus the energy difference becomes 

Ef-Ei = (M z -M z+1 ) + W + q + B(G)-B{A l ) 

= {M z - M z+ i) + W + q+ [B(G) - B(G')] - [B(A') - B(G% (A6) 

In the second line of the equation we have introduced the total electron binding energy for the final neutral atom of 
charge Z in its atomic ground-state configuration, B(G'). 

The Qec value is the difference in the atomic masses of neutral atoms in ground-state configurations: 

Qec = Mz+i(G)-Mz(G') 

= [Mz+i + (Z + l)m e - B(G)] - [M z + Zm e - B(G')] 

= (M z+ i-M z )+m e + [B(G')-B{G)}; (A7) 

and the quantity Wo in Eq. (|A1[) is related to Qec by the equation Wq — Qec ~ m e- Thus, Ef — Ei in Eq. (|A6[) can 
be written as 

E f - Ei = q + W -W + [B(G') - B(A')]. (A8) 
For the energy-conserving delta function we now make a Taylor series expansion about the value q + W — Wq: 

5{E f - Ei) = S(q + W- W ) + S'(q + W- Wo)[B{G') - B(A')} + ... (A9) 



If the first term in this expansion is inserted into the double integral, Eq. (|A3j) . then the expression for / reduces to 
the original form Eq. {AT} since the atomic overlap factor is unity under the assumption that the sum over electronic 
configurations A' can be completed by closure: i.e. J2a> K^'l^)! 2 = J2a'( G \ A ')( A '\ G ) = ( G \ G ) = L The second 
term in Eq. (|A9[) involves a derivative of a delta function. This is handled by an integration by parts, in which the 
rest of the integrand is differentiated with respect to q. No boundary terms survive as the integrand vanishes at the 
boundaries. Thus the atomic overlap correction becomes 

r(Z,W) = l- w ^ w Y.\^ A '\ G )\ 2 W G ')- B ^ 



= l-- f ^ r ^^(G')-X;i(A'|G)| 2 i?(A')J. (A10) 
Next, it is useful to recall the eigenvalue equations satisfied by the atomic states \G) and \A'): 

Z+l Z+l Z+l ^ 

Hi\G) = -B(G)\G) with Hi = ^U-{Z + l)e 2 Y^ - + e 2 ~ > ( AU ) 

i ■ -I 'i r ^ • i " a 

2—1 l — l i<Cj — l J 

Z+l Z+l Z+l 

itj 1 \A l ) = -B{A f )\A f ) with HJ 1 - ~ Ze2 _ + e2 — ' ( A12 ) 

i— 1 i— 1 % i<j—l ■* 

where ti is the kinetic energy of electron i, ri is its distance from the nucleus and r\j is the separation of electrons i 
and j. Note, in particular, that 

z+i „ 

1=1 

Inserting these Hamiltonian expressions into Eq. (|A10[) we obtain 



t ^ n dz 
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W - 
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2 




W - 


W 


2 




W Q ~ 


w 


2 





r(Z,W) = l--—-[B{G') + Y f { G \Hj 1 \A'){A'\G) 



1 - — — \B(G') - (G\Hi — HJ \G) + (G\Hi\G) 

^ l -w^{ B{G ' ) + M B{G) - B{G ^ (A14) 
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TABLE X: Comparison of statistical rate functions calculated without the atomic overlap correction, /without, those calculated 
with it included, /with. The change in the Qec value that would lead to the same change in / is given in the last column. 



Parent 



./w 



/wit 



df/f{%) 



dQ/Q{%) 



dQ (eV) 



io c 

14 Q 

18 Ne 
Mg 
26 Si 

Ar 
Ca 
Ti 



22 



34 
38 
412 



2.30089 
42.7779 
134.484 
418.423 
1029.52 
1967.05 
3414.68 
5327.57 
7040.63 



2.30039 
42.7724 
134.469 
418.386 
1029.44 
1966.91 
3414.46 
5327.24 
7040.21 



0.02178 
0.01277 
0.01093 
0.00877 
0.00767 
0.00707 
0.00647 
0.00612 
0.00597 



0.00436 
0.00255 
0.00219 
0.00175 
0.00153 
0.00141 
0.00129 
0.00122 
0.00119 



83 
72 
74 
72 
74 
77 
78 
81 
84 



T z = 0: 

26m 

34 C1 

38m j^- 

42 Sc 

46y 

50 Mn 
54 Co 
2 Ga 
3 As 
3 Br 
4 Rb 



70 I 
74, 



478.279 
1996.10 
3298.10 
4472.52 
7209.90 
10746.6 
15767.5 
26401.6 
32127.0 
38602.2 
47296.9 



478.237 
1995.96 
3297.88 
4472.24 
7209.47 
10746.0 
15766.6 
26400.2 
32125.3 
38600.1 
47294.5 



0.00880 
0.00711 
0.00663 
0.00643 
0.00598 
0.00565 
0.00537 
0.00557 
0.00545 
0.00539 
0.00423 



0.00176 
0.00142 
0.00133 
0.00129 
0.00120 
0.00113 
0.00107 
0.00111 
0.00109 
0.00108 
0.00105 



75 
78 
80 
83 
84 
86 
89 
102 
104 
107 
109 



where in the last line the order of integration and differentiation has been reversed on the assumption that the binding 
energy as a function of Z behaves in a smooth way. Now, B(G r ) is the electronic binding energy of a neutral atom 
with Z electrons, while B(G) is the same quantity for an atom with Z + 1 electrons. Treating B(G') as a function of 
Z, we can expand B(G') in a Taylor series about B(G): 

B{G') = B(G) - ^B(G) + \^B{G) ... (A15) 
Then, substituting this expression in Eq. (|A14[) , we obtain our final expression for the atomic overlap correction: 

^ = 1 -w^wm B ^- (A16) 



This expression was first obtained by Bahcall [214 

It remains for us to estimate the second derivative of the electronic binding energy of neu tral a toms in their ground- 
state configuration. For this we use binding-energy values from the tables of Carlson et al. [215l |. which were obtained 
from self-consistent Hartree-Fock calculations and have been demonstrated to agree with experimental values to within 
5%. We performed a fit to these tabulated values using a fitting function, aZ b , in three ranges of Z values, with the 
following results: 

13.080Zf 42 eV, 6<Zi< 10 
B(G) = 14.945Z2- 37 eV, 11 < Z, < 30 (A17) 
11.435Z2- 45 eV, 31 < Z, < 39, 
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where Zi is the charge of the parent atom in the beta-decay process. It is conventional to use Z as the charge of the 
daughter nucleus in beta decay; thus for positron decay Zi = Z + 1. The second derivative is easily obtained from 
these expressions. 

We have re-computed the statistical rate function /, with the results being listed in Table 1X1 Those results obtained 
without the atomic overlap correction, Eq. (|A1|) . are given under the heading /without, while those with the correction, 
Eqs. (IA2|) and (|A16|) . are labelled /with- The latter results also appear in column 2 of Table ITXl The fractional 
difference between / w ith and /without in percent is given in column 4 and is of order 0.01%, decreasing with increasing 
mass value. This is a very small correction. Furthermore, the statistical rate function depends on the Q-value to 
the fifth power, so the fractional change in Q that would lead to a change in / of the same size as that induced by 
the atomic overlap correction is even smaller: 1/5 x df //. This percentage change is given in column 5 of Table IXl 
As small as this effect is, it can be seen from the last column of the table that the equivalent change in Q-value 
ranges from 70 to 110 eV, an amount that is similar to the experimental uncertainties on the most precisely measured 
Q- values. 



APPENDIX B: ACCOUNTING FOR SYSTEMATIC ERRORS 
1. Isospin-symmetry- breaking correction, Sc2 

In the past (sec, for example, [f|), we have added a systematic uncertainty to the average corrected Tt value 
to account for an apparent systematic difference between calculations of the isospin-symmetry-breaking correction, 
8c2, that used, on the one hand, Saxon- Woods eigenf unctions a nd, on the other hand, Hartree-Fock eigenfunctions. 
The former method, which was the one used by us [7 , 



method, which was used by Ormand and Brown [181 



178, 



182 



179 



gave consistently larger corrections than the latter 
183(. This was deemed to be a systematic effect caused 



by the different shapes of the Saxon- Woods and Hartree-Fock mean-field potentials. In this section, we re-address 
this issue for two reasons: First, the Hartree-Fock calculations of Ormand and Brown are now 14-24 years old and 
are not available for all the nuclei under study here; and second, our most recent Saxon- Woods calculations Q were 
performed in larger shell-model spaces not matched by the Hartree-Fock calculations. For a valid comparison both 
sets of calculations should be done in identical shell-model spaces. 

The isospin-symmetry-breaking correction, Sc2, for a superallowed [3 transition between + , T=l analog states in 
nuclei with A nucleons is computed from the formula Q 



5c2^ £ s<n<-± J2 s>n>. (bi) 



7r^ , a 



Here S a is the spectroscopic factor for the pickup of a nucleon in quantum state a from an A-particle state of spin 
+ and isospin 1, to an (A — l)-particle state of spin a and isospin T % . There is a complete-set sum over all the 
(A — l)-particle states (called parent states, and denoted ir) in Eq. (jBl|) . The sum is divided into two parts: the 
first is over states with isospin T v = 1/2 and is denoted by tt < ; the second is over states with T v = 3/2, denoted 
7r > . Further, the are radial-mismatch factors, which depend on the difference between the radial wave function 
of a proton bound in the decaying nucleus, u 7 ^ a (r), and that of a neutron bound in the daughter nucleus, u 7 ^ a (r). 
Specifically, the radial-mismatch factors are given by: 



K = / <, a M[< a M-< a (r)] dr 
Jo 

poo 

= 1- < (r)< a (r) dr. (B2) 



Jo 

The radial functions are normalized to J \u(r)\ 2 dr = 1. They are labelled by the parent state tt because their 
asymptotic forms are matched to their separation energies, which in turn depend on the parent state. For example, 
if the parent state is the ground state of the (A — l)-system, then the proton separation energy would be S p and the 
neutron separation energy S n , two quantities given in terms of atomic masses and found in any atomic mass table. If, 
however, tt represents an excited state of the (^4 — l)-system, then the proton and neutron separation energies would 
be S p + E x and S n + E x respectively, where E x is the excitation energy of that parent state. 

To compute &ci from Eq. (jBl|) one needs a set of spectroscopic factors, S a , and a set of radial-mismatch integrals, 
f2 a . The difference between Saxon- Woods and Hartree-Fock calculations lies in the method used to evaluate the 
radial-overlap integrals. In one case, the radial functions u(r) are taken to be eigenfunctions of a Saxon- Woods 
potential; in the other case they are eigenfunctions of a Hartree-Fock mean-field potential. However, for a valid 
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comparison between them the Sc2 must be calculated for both with the same spectroscopic factors, S a . This was not 
done in the past, it being assumed that the model spaces were sufficiently comparable that this would not lead to any 
serious error. However, our most recent Saxon- Woods calculation [3] have significantly increased the model spaces 
over the ones used before, so now it is essential that a new Hartree-Fock calculation be undertaken. We report such 
a calculation here. 

To be clear about the procedure, we illustrate it for the specific case of the decay of 34 C1 to 34 S. The decaying 
nucleus has Z + 1 = 17 protons; the daughter nucleus has Z=16 protons. In the Saxon- Woods approach, the proton 
radial wave functions are taken to be eigenfunctions of a potential defined for a nucleus of mass A and charge Z + 1 
as follows: 



V(r) = -V f(r) - V 8 g{r)\.a + V c (r), 



(B3) 



where 



f(r) = {1 + exp ((r-i?)/a)} \ 
9{r) = 



m-nC J a«r 



x < 1 + cxp 



1 (r-R> 
exp 



a s 
r - R t 
a.. 



V c (r) = Ze 2 /r, for r > R c 
Ze 2 ( r 2 \ 



(B4) 



Here, R = r (A - l) 1 / 3 and R s =r s (A- l) 1 / 3 . The three terms in Eq. (|B3j are the central, spin-orbit and Coulomb 
terms respectively. In our calculations Q most of the parameters were fixed at standard values, with the well depth Vq 
being adjusted case by case so that the binding energy of the eigenfunction being computed matched the separation 
energy to the corresponding parent state - in 33 S, for our example. Likewise the neutron radial functions were taken 
to be eigenfunctions of a similar potential but with the Coulomb term omitted. 

Although the Hartree-Fock procedure is comparable, there is one issue unique to this approach which requires 
particular attention. For our illustrative example, a Hartree-Fock calculation might first be mounted for 3 CI, which 
would yield a mean field with central, spin-orbit and Coulomb terms. The required proton radial functions would 
then be taken as eigenfunctions of this mean field with the strength of the central term readjusted case by case so 
that the computed binding energy matched the appropriate separation energy. A second Hartree-Fock calculation 
might then be mounted for 34 S, from which the neutron radial functions would be similarly determined in the mean 
field, but without the Coulomb term. However, under these circumstances, if the Coulomb terms in the Hartree-Fock 
mean-field potential were to be compared with those in the Saxon- Woods potential, a very significant difference would 
emerge. In the Hartree-Fock case, the Coulomb term is 



V c (r) 



d 3 r 



r — r 



3e*_ 
2 



:P P (i 



1/3 



(B5) 



which depends on the proton density (of 34 C1 in our example) that is generated as part of the Hartree-Fock procedure. 
The two terms in Eq. (IB5|) are called the direct and exchange terms respectively. If we take the asymptotic limit of 
the direct term for large r, we obtain 



r) 



'3„' 



r — r 



e 
r 



d 3 rV P (r') 



(Z + l)e 2 



(B6) 



Since the Hartree-Fock proton density is normalized to (Z + 1) protons in CI, the asymptotic form of the Coulomb 
potential tends to (Z + l)e 2 /r. However, this disagrees with the equivalent Saxon- Woods calculation, which has the 
form Ze 2 /r (see Eq. lB4|) . 

This discrepancy is important and constitutes, in our opinion, a serious flaw in this Hartree-Fock calculation of 
the radial-mismatch factor. Since a proton removed from a nucleus of charge Z + 1 leaves behind Z protons, its 
asymptotic interaction is with charge Z - as described by the Saxon- Woods potential - and not with charge Z + 1 . 
This deficiency in Hartree-Fock would be cured in principle by the Coulomb exchange term. However, in Skyrme- 
Hartree-Fock calculations it is not possible to compute the exchange term exactly without sacrificing the simplicities 
that come with use of zero-range Skyrme interactions. The exchange term appearing in Eq. (|B5|) is a commonly 
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TABLE XI: Adopted isospin-symmetry-breaking corrections, 8c2 in percent units, and their assigned uncertainties obtained 


from Hartree-Fock calculations. Also listed are earlier results obtained with Saxon- Woods (SW) eigenfunctions, as published 


m [7|J. 






Nucleus 


HF 


SW 


T z = -1 : 






10 C 


0.215(35) 


0.165(15) 


14 Q 


0.255(30) 


0.275(15) 


18 Ne 


0.205(55) 


0.410(25) 


22 Mg 


0.250(55) 


0.370(20) 


Si 


0.335(55) 


0.405(25) 


30 q 


11.0411(00 J 


U. 1 1)11(211 ) 


34 a 

Ar 


O.OlU(OU) 


0.000(00) 


L>a 


U.DULUbU ) 


U. I 40( ( U ) 


42 rp- 
1 1 


U.0o0(0U J 


U.oo0( / OJ 


T z = : 






26 Al 


0.410(50) 


0.280(15) 


34 C1 


0.595(55) 


0.550(45) 


38 R 


0.640(60) 


0.550(55) 


42 Sc 


0.620(55) 


0.645(55) 


46y 


0.525(55) 


0.545(55) 


50 Mn 


0.575(55) 


0.610(50) 


54 Co 


0.635(55) 


0.720(60) 


62 Ga 


0.93(16) 


1.20(20) 


66 As 


1.10(35) 


1.35(40) 


70 Br 


1.14(25) 


1.25(25) 


74 Rb 


1.29(16) 


1.50(30) 


used local approximation 


, which might well be appropriate for the nuclear interior and for the computation of bulk 


properties such as bindinj 


5 energies and radii, but it certainly does not do the job asymptotically, which is 


the region 



of greatest importance to our calculations. 

To circumvent this difficulty, we have chosen to alter the Hartree-Fock protocol. Instead of mounting two Hartree- 
Fock calculations - for 34 C1 and 34 S - as just described, we mount a single calculation for the nucleus with (A — 1) 
nucleons and Z protons - 33 S in our example. We then use the proton mean held from this calculation to generate 
the proton eigenfunctions, a (r); and the neutron mean held from the same calculation to generate the neutron 
eigenfunctions u^ a (r). In this procedure, the Coulomb interaction automatically has the correct asymptotic form. 
It is also fully consistent with the Saxon- Woods potential parameterisation, Eq. (|B4|) . which considers the nucleus 
of mass (A — 1) as the core to which the last particle is bound, since the radius of the potential is parameterised as 
Tq(A — l) 1 / 3 rather than tqA 1 / 3 . Calculations of Sc2 with this new Hartree-Fock protocol will be presented in the 
next section. It can be noted here, however, that these results are larger than those obtained with the conventional 
protocol by between 10% to 40% depending on the Skyrme interaction used and the nucleus under study. This change 
of protocol goes a long way in reducing the systematic error between Saxon- Woods and Hartree-Fock calculations. 

2. New Hartree-Fock calculations for 5c2 

Here we present our new Hartree-Fock calculations for the isospin-symmetry-breaking correction, 5c2, for the 20 
cases of superallowed Fermi beta decay considered in our survey: they range from 10 C to 74 Rb. Our procedure was 
that discussed at the end of the last section and involved obtaining the mean held from a Hartree-Fock calculation 
in the (A — l)-system. The proton and neutron radial functions were obtained as eigenfunctions of this mean-held 
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potential, whose overall strength was scaled on a case-by-case basis to ensure the eigenfunction's asymptotic solution 
matched the required separation energy These cigenfunctions were used to compute the radial-mismatch factors, f2^ 
of Eq. (|B2jl . For the spectroscopic factors needed in Eq. (|B1|) , we ran several shell-model calculations with the model 
spaces and effective interactions used recently M] in calcula tions with Sax on- W oods potentials. We also considered 
three choices of the Skyrme interaction: S GII 2161 . SkM* |217l ] and Ska [218]. The first two of these interactions 
were the ones used by Ormand and Brown [182 . Il83| in their computations of 8c2- The third interaction, Ska, is of 
similar quality and was one of the first to be fitted to the incompressibility of nuclear matter, a key constraint used in 
all later Skyrme interactions. More re cent Skyrme in teractions tend to be used in conjunction with pairing forces in 
Hartree-Fock-Boglioubov calculations [219l l22dl l22l| . Since we have not included pairing forces in the present work, 
we have not attempted to use any of these more recent interactions. 

The results from our Hartree-Fock calculations for b~c2 are listed in column 2 of Table IXIl For each transition, 
the central value is an average of the results obtained with the three choices of Skyrme interactions. To assign an 
uncertainty, we have examined the spread in results obtained with the different Skyrme interactions and with different 
shell-model effective interactions and model spaces. We also list in the third column of the table the values we adopted 
from our Saxon- Woods computations, which originally appeared in ref. Q. 
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